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ABSTRACT 

We present a comprehensive catalog and analysis of broad-band afterglow observations for 103 short-duration 
gamma-ray bursts (GRBs), comprised of all short GRBs from November 2004 to March 2015 with prompt 
follow-up observations in the X-ray, optical, near-infrared and/or radio bands. These afterglow observations 
have uncovered 71 X-ray detections, 30 optical/NIR detections, and 4 radio detections. Employing the standard 
afterglow synchrotron model, we perform joint probability analyses for a subset of 38 short GRBs with well- 
sampled light curves to infer the burst isotropic-equivalent energies and circumburst densities. Eor this subset, 
we find median isotropic-equivalent 7 -ray and kinetic energies of £ 7 ,iso ~ 2 x 10^' erg, and Ek, iso ~ (1 - 3) x 
10^* erg, respectively, depending on the values of the model input parameters. We further find that short GRBs 
occur in low-density environments, with a median density of « « (3- 15) x 10“^ cm“^, and that « 80-95% 
of bursts have densities of n < 1 cm“^. We investigate trends between the circumburst densities and host 
galaxy properties, and find that events located at large projected offsets of > 10 effective radii from their hosts 
exhibit particularly low densities of n < 10“^ cm“^, consistent with an IGM-like environment. Using late-time 
afterglow data for 11 events, we find a median jet opening angle of 9j = 16 ± 10°. We also calculate a median 
beaming factor of ft, « 0.04, leading to a beaming-corrected total energy release of Etme « 1.6 x 10^*^ erg. 
Eurthermore, we calculate a beaming-corrected event rate of JRtme = 270^}gg° Gpc“^ yr“', or « 8 ^ 5 ^ yr“^ within 
a 200 Mpc volume, the Advanced LIGO/Virgo typical detection distance for NS-NS binaries. 


1. INTRODUCTION 

The afterglow emission from gamma-ray bursts (GRBs) 
provides a unique way to probe their basic properties: 
the energy scale, circumburst densities, and jet opening 
angles. Nearly two decades of long-duration GRB after¬ 
glow observations established a median beaming-corrected 
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In some cases, the radial density 
profiles of their circumburst environments r eflect the wind 
environments expecte d for massive stars (IChevalier & Lil 
[ToOOtlYost etlklU?)?)! . 
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for providing precise burst localiza tion, and therefor e robust 
associations to host galaxies (e.g., lEong et ^l2013h . These 
host studies have revealed that at least some short GRBs orig¬ 
inate from older stellar populations, and are di stinct in their 
global properties from t he hosts of long GRBs (lBergeijl2009l: 
iLeibler & Bergeiil2010l: lEong et al.ll2013l: lBergeill2014l). The 
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the lack of associated supernovae (e.g., lE ox et al. l 120051: 

ISoderberg et al.l2006l:lKocevski et aklioiOtlFong et al.l2014ll . 

and the discov ery of a near-IR “kilonova” from the short 
GRB 130603B (iBerger et al.l 1201 31: iTanvir et al.l IToTl . have 
provided strong observational support for a compact object 
binary merger progenitor involving two neutron stars or a 
neutron star and a stellar mass black h ole (NS-NS/NS-BH; 
lEichler et akllT^lNaravaner^fT^ . 

As studies of the host galaxies of short GRBs have pro¬ 
gressed at a rapid pace, knowledge of their basic explo¬ 
sion properties has been limited by both the paucity of af¬ 
terglow detections and the relatively low detection rate of 
well-localized (^few arcsec uncertainty) short GRBs from 
the Swift satellite. Eurthermore, afterglow studies greatly 
benefit from observations across several orders of frequency, 
which serve to provide significantly tighter constraints on 
the energy scales and circumburst densities than single 
band observations. Thus far, short GRB afterglow stud¬ 
ies have either foc used on a small number of events in 
a single band (e.g., iBergeil l2()07l: iNakad 1200'^ iKann et akl 
l2008t iNicuesa Guelbenzu et al.ll2012a^ , or on radio through 
X-ray observations of individual burst s (IBerger et al.1 l2005t 
ISoderberg et al.ll2()()^lFong et al.ll2014ft . Comparative studies 
relative to long GRBs have only served to argue for lower en¬ 
ergy scales and circumb urst densities , but have not provided 
actual distributions (e.g.,|Berge320T3)- Thus, there have been 
no attempts to utilize the full set of broad-band afterglow ob¬ 
servations of short GRBs. 

The relative difficulty of detecting short GRB afterglows 
is likely a direct reflection of a combination of low explo¬ 
sion energies and low circumburst densities. Predictions 
for NS-NS/NS-BH mergers span several orders of magni¬ 
tude in density, depending on the precise distribution of na- 
tal kick velocities, merger t imescales, and host gal axy type 
(iPerna & Belczyns^ l2002t iBelczynski et al.l l2006h . Simi- 
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larly, different mechanisms of energy extraction to power a 
relavistic blast-wave can produce en ergy scales which differ 
by three order s of magnitude (e.g., [Ruffert & Jankalfl999bl: 
lRosswogl2005h . Furthermore, the calculation of the true short 
GRB energy scale and event rate directly depends on the ge¬ 
ometry of their jets. Constraints on the energetics, densities, 
and jet opening angles from short GRB afterglows thus of¬ 
fer a way to study these fundamental questions. These ba¬ 
sic properties also serve as critical inputs for the detectabil¬ 
ity of other electromagnetic counterparts to compact object 
mergers, and will directly affect follow-up strategies to gravi- 
tatio nal wave signals detected with Advanced LIGO /VIRGO 
(e.g.. lMetzger & Bergeiil201^lNissanke et al.ll2013h . 

Taking advantage of the dedicated searches for short GRB 
afterglows at all wavelengths, we are now in a position to ex¬ 
plore these basic properties for a large population of events 
for the first time. Here, we present the first comprehensive 
broad-band catalog of short GRB afterglows, representing a 
decade of observations since the launch of Swift in 2004, and 
use this sample to constrain short GRB energies and circum- 
burst densities. In Section |2] we introduce the sample and 
data reduction methods for X-ray through radio observations. 
In Section [3 we model the temporal and spectral behavior of 
the afterglows, and use the observations to infer constraints 
on the energetics and circumburst densities for each burst. In 
Section |4l we present the energetics and circumburst densi¬ 
ties for 38 events with well-sampled afterglow data sets. In 
Section 13 we discuss the observational afterglow properties, 
jet opening angles, and investigate trends between the bursts 
and their galactic environments. Finally, in the Appendix, we 
provide X-ray, optical, and radio afterglow catalogs. In an up¬ 
coming work, we will address the effect of the basic inferred 
properties on the detectability of electromagnetic counterparts 
to gravitational wave sources. 

Throughout the paper, all magnitudes are in the AB sys¬ 
tem and are corr ected for Galactic extinction in the direc- 
tion of each burst (ISchlegel et al.ll9^ISchlaflv & Finkbeineii 
1201 Ih . Unless otherwise noted, uncertainties correspond to 
Icr confidence. We employ a standard ACDM cosmology 
with Qm = 0.27, = 0.73, and T/q = 71 km s“' Mpc“'. 

2. OBSERVATIONS 
2.1. Sample 

We present afterglow observat ions for 103 short GR Bs dis¬ 
covered by the Swift satellit e (iGeh rels et a ll l2004ll . Fermi 
satelli te (lAtwood et akl l2009t iMeegan et al.l l2009l) . Konus- 
Wind (I Aptekar et al.ll995lk High Energy Transient Explorer 2 
(HET E-2; [Ricker et alJl200l . or the Interplanetary Network 
(IPN: [Hurley et al.l2010h between November 2004 and March 
2015. We restrict our sample to all bursts with Tgo <2 s 
and follow-up observations in any of the X-ray, optical, near- 
infrared (NIR) or radio bands on timescales of 6t < few days 
(where 6t corresponds to the time after the 7-ray trigger). 
We also include three events (GRBs 050724A, 090607 and 
100213A) which have Tgo ~ 2.5-3 s but which exhibit the 
spectra l hardness and neglig i ble spectral lags ty pical of short 
GRBs dMarshall et al.ll200^ iGrupe et al.ll20f(^ . Eor bursts 
with optical/NIR follow-up, we only include data from bursts 
with afterglow detections or meaningful limits of > 20 mag at 

< 1 day. Basic information for the sample of 103 events, 
including durations, redshifts, and the available follow-up in 
each observing band is presented in Table [T] 

Of the 103 bursts in our sample, 96 (93%) have X-ray 


Table 1 

Short GRB Basic Information 


GRB 

Tgo 

Z 

X-ray 

Opt/NIR 

Radio 


(s) 





050202 

0.27 



N 

N 

050509B 

0.04 

0.225 

Y 

N 

N 

050709 

0.07 

0.161 

Y 

Y 

N 

050724A 

3.0 

0.257 

Y 

Y 

Y 

050813 

0.6 

0.72/1.8 

Y 

N 

N 

050906 

0.26 

N 

N 

N 

050925 

0.07 


N 


N 

051105A 

0.09 


N* 


N 

051210 

1.3 

> 1.4 

Y 

N 


051221A 

1.4 

0.546 

Y 

Y 

Y 

060121 

2.0 

<4.1 

Y 

Y 


060313 

0.7 

< 1.7 

Y 

Y 

N 

060502B 

0.09 

0.287 

Y 

N 


060801 

0.5 

1.130 

Y 

N 

N 

061006 

0.4 

0.438 

Y 

Y 


061201 

0.8 

0.111 

Y 

Y 


061210 

0.2 

0.41 

Y 

N 

N 

061217 

0.2 

0.827 

Y 

N 


070209 

0.09 


N 

N 


070406“ 

1.20 


N 

N 


070429B 

0.5 

0.902 

Y 

N 

N 

070707“ 

1.1 

<3.6 

Y 

Y 


070714B 

2.0 

0.923 

Y 

Y 

N 

070724A 

0.4 

0.457 

Y 

Y 

N 

070729 

0.9 

0.8 

Y 

N 

N 

070809 

1.3 

0.473 

Y 

Y 


07081OB 

0.08 


N 

N 


070923 

0.05 




N 

071017“ 

0.5 


N 



071112B“ 

0.30 


N 

N 

N 

071227 

1.8 

0.381 

Y 

Y 


080121" 

0.7 


N 

N 


080123 

0.4 


Y 

N 


080426 

1.7 


Y 

N 


080503 

0.3 

<4.2 

Y 

Y 

N 

080702A 

0.5 


Y 

N 

N 

080905A 

1.0 

0.122 

Y 

Y 


080919 

0.6 


Y 

N 


081024A 

1.8 


Y 

N 

N 

081024B" 

0.4 


N* 

N 

N 

081226A 

0.4 

<4.1 

Y 

Y 

N 

081226B" 

0.7 


N 

N 

N 

090305A 

0.4 

<4.1 

Y 

Y 


090417A 

0.07 




N 

090426 

1.3 

2.609 

Y 

Y 


090510 

0.3 

0.903 

Y 

Y 

N 

090515 

0.04 

0.403 

Y 

Y 

N 

090607 

2.3“ 


Y 

N 


09062IB 

0.14 


Y 

N 

N 

090715A 

0.5 



N 

N 

090916“ 

0.3 


N 



091109B 

0.3 

<4.4 

Y 

Y 


091117“ 

0.43 


N* 

N 

N 

100117A 

0.30 

0.915 

Y 

Y 


100206A 

0.1 

0.407 

Y 

N 


100213A 

2.4'' 


Y 



100625A 

0.3 

0.452 

Y 

N 

N 

100628A 

0.04 


N 

N 

N 

100702A 

0.16 


Y 

N 


101219A 

0.6 

0.718 

Y 

N 


101224A 

0.2 


Y 


N 

110112A 

0.5 

<5.3 

Y 

Y 

N 

110112B" 

0.5 


N 

N 

N 
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Table 1 — Continued 


GRB 

Tgo 

(s) 


X-ray 

Opt/NIR 

Radio 

110420B‘' 

0.08 


N* 

N 

N 

111020A 

0.4 


Y 

N 

N 

111117A 

0.5 

1.3 

Y 

N 

N 

111121A 

0.45 


Y 


N 

111222A“ 

0.3 


Y 



120229A 

0.22 



N 

N 

120305A 

0.1 


Y 

N 

N 

120521A 

0.45 


Y 

N 

N 

120630A'' 

0.6 


Y 



120804A 

0.81 

1.3 

Y 

Y 

N 

120817B“ 

0.19 


N 

N 


121226A 

1.0 


Y 

N 

N 

130313A 

0.26 


Y 

N 

N 

130515A 

0.29 


Y 

N 


130603B 

0.18 

0.356 

Y 

Y 

Y 

130626A 

0.16 


N 



130716A 

0.8 


Y 

N 

N 

130822A 

0.04 


Y 

N 

N 

130912A 

0.28 

<4.1 

Y 

Y 

N 

131004A 

1.54 

0.717 

Y 

Y 

N 

131125A^ 

0.5 



N 


131126 A^ 

0.3 



N 


131224A“ 

0.8 


N 


N 

140129B 

1.36 

< 1.5 

Y 

Y 


140320A"' 

0.45 


Y 

N 


140402A“ 

0.03 


N 

N 


140414A“ 

0.7 


N 

N 


140516A 

0.19 


Y 

N 

N 

140606A“ 

0.34 


N 

N 


140619B“ 

0.5 


N 

N 

N 

140622A 

0.13 

0.959 

Y 

N 

N 

140903A 

0.30 

0.351 

Y 

Y 

Y 

140930B 

0.30 

<4.1 

Y 

Y 

N 

141202A“ 

1.3 


N 



141205A" 

1.1 


N 

N 


141212A 

0.30 

0.596 

Y 

N 

N 

150101A 

0.06 


Y 


N 

150101B 

0.018 

0.134 

Y 

Y 

N 

150120A 

1.20 

0.460 

Y 

N 

N 

150301 A" 

0.48 


Y 




Note. — Short GRBs with X-ray, optical/NIR or radio follow¬ 
up observations. “Y” = detection, “N” = non-detection, and 

• • • means there is no follow-up in that band. 

“ Observing constraint or delayed Swiftl'KSH observations. 

* There are XRT observations but no X-ray flux upper limit is 
reported. 

This burst has Tqo = 2.3 ±0.1 s but is spectrally hard 
iMarshall et alJ2009l) . 

This burst has 7% = 2.4 ± 0.4 s, but the spect ral lag of 5 ± 15 ms 
indicates this is a short-hard burst dGrune et alJl201(]h . 

Delayed reporting of burst position preventing immediate ground- 
based follow-up. 

^ IPN-localized burst with no Swift follow-up. 

follow-up observations, 87 (84%) have optical/NIR observa¬ 
tions, and 60 (58%) have radio observations. These obser¬ 
vations have uncovered 71 X-ray afterglows, 30 optical/NIR 
afterglows, and 4 radio afterglows, leading to detection frac¬ 
tions of 74%, 34%, and 7%, respectively (Table [T]). The ob¬ 
servations for these bursts are catalogued in the Appendix (Ta¬ 
bles A1-A3). 

Thirty-one bursts in this sample have redshifts deter¬ 


mined from their host galaxies (Table [T]i, while two 
events have spectroscopic redshift s from absorp ti on in 
their afterglows, G RB 090426A: jAntonell i et alj 1 20091: 
Levesque etal.l I2OIOI) and GRB 13 0603B dCucchiara et alJ 
20131 Ide Ugarte Postigo et alJl2014ft . Furthermore, we place 
upper limits on the redshift from the detection of the optical 
afterglow, and therefore the lack of suppression blueward of 
the Lyman limit (Aq = 912 A) or Lyman-a line (Aq = 1216 A), 
for 11 bursts (Table [U. We also place a lower limit of 
z > 1.4 for GRB 051210 from the lack of detection of emis- 
sion or absorption f eatures in the spectrum of the host galaxy 
(iBerger et al.ll200'^ . In addition to the broad-band afterglow 
observations that have been published thus far, we present 
new optical/NIR observations for 11 bursts (Table A2), and 
new radio observations for 25 events (Table A3). 

2.1.1. Observing Constraints 

Of the 25 bursts with X-ray observations and no detected 
X-ray afterglow, 18 events had a delayed SwiftlXKf response 
due to an observing constraint or burst discovery from an¬ 
other satellite (indicated in Table [T]); thus, the non-detection 
of X-ray afterglows for these events are due to factors unre¬ 
lated to the bursts themselves. Likewise, 12 bursts with op¬ 
tical observations and no detected optical afterglow have an 
observing constraint (e.g., delayed precise localization from 
Swift, crowded field, high Galactic extinction sightline, po¬ 
sition contaminated by a bright star), making the detection of 
an optical afterglow particularly challenging (Table A2). Only 
one burst (GRB 150101B) with radio observations is severely 
contaminated by a neighboring bright source, thus preventing 
a deep radio limit at the afterglow position. Taking these ob¬ 
serving constraints into account, we find that the vast majority 
(91%) of bursts that have X-ray follow-up and no observing 
constraints result in an X-ray afterglow detection. Similarly, 
the fraction of detected optical afterglows increases to 40% 
after correcting for observing constraints. The radio detection 
fraction remains at r; 7%. 

2.2. X-rays 

We gather all available X-ray a fterglow observations from 
the Swift light curve repositorjQ (lEvans et al.l 12007 al l2009h . 
the GRB Coo rdin ates Network (GCN) Circulars, and the lit¬ 
erature (Table IaTI) . The data were taken with the X-ray Tele¬ 
scope (XRT) on-board Swift, the Chandra X-ray Observatory, 
and the X-ray Multi-Mirror Mission (XMM-Newton). Ten 
bursts have Chandra observations , wh ile three bursts have 
XMM-Newton observations (Table IaTTi . We use unabsorbed 
fluxes and uncertainties in the 0.3 -10 keV energy band when 
they are available; otherwise, we use the count-rate light 
curves in the same energy range and convert to fluxes (de¬ 
scribed later in this section). For bursts with multiple upper 
limits, we only include those which help to constrain the tem¬ 
poral behavior of the X-ray light curve. Of the 96 short GRBs 
with X-ray observations, four events have no reported mea¬ 
surements or upper limitfl Therefore, the resulting late -tim e 
X-ray afterglow catalog is comprised of 92 events (Table lATI) . 

When applicable, we convert the count rate light curves 
to unabsorbed fluxes using the count-rate-to-unabsorbed-flux 
conversion factors provided by the Swift light curve reposi¬ 
tory. These factors are derived from the automatic spectral 

^ http://www.swift.ac.uk/xrt_curves 

5 These events are GRBs051105A, 081024B, 091117, and I10420B. 
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Figure 1. Broad-band afterglow light curves of all short GRBs with follow-up observations between November 2004 and March 2015. In all panels, circles 
denote detections, triangles indicate 3 (t upper limits, and solid lines connect observations for the same burst. Top left: 0.3 — 10 keV X-ray afterglow light curves 
for 92 short GRBs with X-ray observations con'esponding to the events in Table Al. Included are 3cr upper limits for 29 bursts with no detected X-ray emission 
at 5t > 1000 sec (blue triangles). The data for 5t > 1000 sec (grey vertical band) are used in our subsequent afterglow analysis. Top right: Optical afterglow 
observations for 87 short GRBs corresponding to the events in Table A2. The 30 bursts with detected optical afterglows are shown (orange circles), along with 
3cr upper limits for bursts with an X-ray afterglow detection (hlled triangles), and bursts with no X-ray afterglow detection (open triangles). Bottom left: Radio 
afterglow data for 60 short GRBs corresponding to the events in Table A3. The light curves for four short GRBs with radio afterglow detections are shown: 
GRBs 050724A (black), 051221A (red), 130603B (light red), and 140903A (maroon). Upper limits for the remaining 56 events with no detected radio afterglows 
are shown (3cr; grey triangles), including 25 bursts for which radio observations have not been published in the literature thus far. 


fitting routine (lEvans et al.ll2009t) . This routine fits the X-ray 
spectrum for each burst to an absorbed power law model char¬ 
acterized by photon index, F, and the intrinsic neutral hydro¬ 
gen absorption column, A^H.int, in excess of the Galactic col¬ 
umn density in the direction of the burst (iKalberla et all2005l : 
IWakker et al.ll20lT| ; IWillingale et alj|2013h . We use spectral 
parameters extracted in the Photon Counting (PC) mode when 
available; otherwise, we use parameters from the Windowed 
Timing (WT) mode. In ten cases, the value of A^n.int is highly 
uncertain, but consistent with zero. Therefore, utilizing the 
median value of A^H.int niay result in an over-estimate of the 
true unabsorbed flux. Instead of using the given conversion 
factors for these bursts, we calculate the unabsorbed fluxes us¬ 
ing WebPIMMlfl setting A^H.int = 0. For 16 events, no count- 
rate-to-unabsorbed-flux conversion factor is available, so we 

^ https://heasarc.gsfc.nasa.gov/Tools/wSpimms.html 


employ a fiducial value of 1 x 10“" erg cm“^ s“' cf' set by 
the median value for all of the events in our sample. Applying 
these conversion factors to each of the count-rate light curves 
from SwiftlXKY, Chandra and XMM-Newton, we obtain the 
unabsorbed fluxes, Icr uncertainties, and 3 (t upper limits for 
each burst (Table Al). 

To enable comparison of the X-ray light curves to the op¬ 
tical and radio data, we convert the X-ray fluxes to flux den¬ 
sities, F^ x, at a fiducial energy of 1 keV (/v,x oc where 
/3x = 1 -F). When no spectral information is available, we 
use a fiducial spectral index of px.med = -TO, set by the me¬ 
dian value of the events in our sample. The flux densities, Icr 
uncertainties, and 3cr upper limits are listed in Table Al, and 
the resulting light curves are shown in Figure [T] 

Finally, we present Chandra observations for three bursts 
that have not been published in the literature thus far: 
GRBsl20804A (PI: Burrows), 140930B (PI: Fong) and 




















Short GRB Afterglows 


5 


150101B (Pis: Trqja, Levan). We retrieve the pre-processed 
Level 2 data from the Chandra archive. We use the Cl AO 
data reduction package to extract a count-rate within a 2.5"- 
radius source aperture centered on the X-ray afterglow posi¬ 
tion, and utilize source-free regions on the same chip to esti¬ 
mate the background. For GRB 120804A, we obtain spectral 
parameters from ea rlier Chandra epochs of the same burst 
(lBergeretal.ll2^ . For GRBs 140930B and 150101B, we 
use CIAO/specextract to extract a spectrum and obtain 
the spectral parameters. We then apply these parameters to 
the count rate to convert to flux density as described above 
(Table Al). 

2.3. Optical/NIR 

For each burst we gather all available optical and NIR af¬ 
terglow observations from the literature and GCN Circulars 
(see Table A2 for references). When there are multiple up¬ 
per limits for a given burst, we include only the deepest lim¬ 
its at (5f < 1 day. We convert all magnitudes to the AB sys¬ 
tem using instrument-specific conversion factors when avail¬ 
able, o r the standard conversions following iBlanton & RoweisI 
(120071) . We correct all observations for Galactic extinc¬ 
tion in the direction of ea ch burst (iSchlegel et alJ 119981 
iSchlaflv & Finkbeineill201 ih . and convert AB magnitudes to 
flux densities, /v.opt- A log of observations for 87 events with 
optical/NIR follow-up is provided in Table A2, and the light 
curves and upper limits are shown in Figure[T] Bursts with no 
detected optical afterglow are further classified by the detec¬ 
tion of an X-ray afterglow (Figure [1]). 

We also present optical/NIR observations of 11 short 
GRBs that have not been published in the literature thus far: 
GRBs070724A, 100628A, 110420B, 120229A, 130716A, 
140402A, 140606A, 140619B, 140930B, 150101B, and 
150120A. These observations were enabled by target-of- 
opportunity programs on the twin 6.5-m Magellan tele¬ 
scopes (PI: Berger), the twin 8-m Gemini telescopes (Pis: 
Berger, Cucchiara, Fox, Tanvir), and the 6.5-m MMT (Pis: 
Berger, Fong). We use standard tasks in the IRAF/ccdred 
package to process the Magellan and MMT data, and the 
IRAF/gemini package to process the Gemini data. For each 
of these bursts, we obtained the first epoch of observations at 
(5f Ri 1 - 20 hr and at least one additional set of observations 
at St > 24 hr to provide a template (Table A2). To assess any 
fading between the two epochs, we perform digital ima ge sub¬ 
tracti on for each burst and filter using the ISIS software (lAlardI 
[2000h . With the exception of GRBs 070724A, 140930B, and 
150101B, we find no significant emission in any of the sub¬ 
tracted images. We therefore employ aperture photometry us¬ 
ing standard tasks in IRAF to place 3cr upper limits on the 
optical/NIR afterglow brightness. To measure the afterglow 
brightness for GRBs 070724A and 140903B, we employ aper¬ 
ture photometry on the detected point source in the subtracted 
image. For GRB 150101B, the afterglow position in the sub¬ 
tracted image is contaminated by residual emission from its 
bright host galaxy; thus, we use PSF photometry to measure 
the afterglow brightness. Details of the photometry will be 
outlined in an upcoming paper (Fong et al., in prep.) The ob¬ 
servational details, afterglow brightness and 3cr upper limits 
for these bursts are presented in Figure [1] and listed in Ta¬ 
ble A2. 

2.4. Radio 

We gather all available radio afterglow data taken with the 
Karl G. Jansky Very Large Array (VLA), Westerbork Syn¬ 


thesis Radio Telescope (WSRT), Australia Telescope Com¬ 
pact Array (ATCA), and Combined Array for Research in 
Millimeter-wave Astronomy (CARMA). The resulting radio 
afterglow catalog is comprised of 60 short GRBs (Table A3). 
The large majority of events, 53 (88%), were observed with 
the VLA (Table A3), 28 of which were observed with the up- 
graded VLA, whic h has a ten-fold improvement in sensitivity 
(IPerlev et al.ll201 it) . We present new observations enabled by 
target-of-opportunity programs with the upgraded VLA (PI: 
Berger) and CARMA (PI: Zauderer) for 25 bursts (Table A3). 
In 11 cases, we obtained multiple sets of observations to probe 
the radio emission on timescales spanning (5f Ri 1 - 10 days. 

For data calibration and analysis of VLA observations, we 
follow standard procedu res in the Astr onomical Image Pro¬ 
cessing System (AIPS: iGreisenI [200^ . For CARMA ob¬ 
servations, we use the Multichannel Image Reconstruction, 
Image Analysis a nd Display (MIRIAD) software package 
(iSault et alJ Il995h . For the majority of cases, we do not 
find any uncatalogued radio sources in or around the X- 
ray or g-ray positions. To calculate the 3cr upper limits on 
the radio afterglow brightness, we measured the RMS noise 
in the map using a large source-free central region utiliz¬ 
ing AIPS/IMSTAT for VLA data and MIRIAD/IMSTAT for 
CARMA data. The radio afterglow detections and upper lim¬ 
its for 60 short GRBs with radio observations are listed in 
Table A3 and displayed in Figure [T] Flux densities reported 
here supercede those reported in GCN Circulars. 

Four bursts have detected radio afterglows, an d all discov¬ 
eries w ere made w i th the VLA (GRB 05072 4A:lBerger et al. 
2001 iPanaitescul [200l GRB 05122 1 A: iSoderberg et~ 
Boa GRB 130603B nBng et akllml GRB 140903A: this 
work). In particular, we present the discovery of the radio af¬ 
terglow of GRB 140903A, which is detected at three frequen¬ 
cies: 1.4, 6.0 and 9.8 GHz (Figure [1] and Table A3), and will 
be further discussed in an upcoming work. 

2.5. Afterglow Brightness and TOO Response Times 

Thanks to Swift, X-ray and optical follow-up typically com¬ 
mences within Ri 60 sec after the GRB is detected. While 
SwiftlXKY is responsible for nearly all of the X-ray afterglow 
detections, only three short GRBs have detected optical af¬ 
terglows with 5w(/f/UVOT (Table A2). Thus, the detection 
of short GRB afterglows, or placement of meaningful upper 
limits, in the optical and radio bands relies on ground-based 
Target-of-Opportunity programs. For the 30 short GRBs 
with detected optical afterglows, the median optical afterglow 
brightness is Ri 23.0 mag at St Ri 7.0 hr. The median 3cr limit 
placed on bursts with an X-ray afterglow and no detected op¬ 
tical afterglow is >23.8 mag with a median response time of 
St Ri 7.4 hr, while the limit placed on bursts with no detected 
X-ray or optical afterglow is less constraining and more de¬ 
layed, > 22.7 mag at St Ri 12.2 hr (Figure [TJ. This is likely 
due to the more limited number of facilities with instruments 
that can cover the comparatively larger g-ray positional error 
circles. However, since we only include bursts with optical 
limits of > 20 mag, we are excluding a fraction of bursts with 
very shallow follow-up; thus the limits here for bursts with 
only g-ray localizations are an optimistic representation of 
the entire population. In the radio band, the median 3 (t up¬ 
per limit for all observations is < 74.1 ply, and the median 
response time to the first observation is St k, 24.7 hr. In more 
recent cases, we set unprecedented 3 (t limits of 15-20p,Jy 
using the upgraded VLA on timescales of (5f Ri 1 - 10 days. 
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3. BROAD-BAND AFTERGLOW ANALYSIS 
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We utilize the broad-band afterglow observations to con¬ 
strain the explosion properties and circumburst environment 
of each burst. We adopt the standard synchrotron model for a 
relativistic blastwave in a c onstant density medium (ISari et alJ 
ll998t[Granot & Saril2002l) . as expected for a non-massive star 
progenitor. This model provides a mapping from the broad¬ 
band afterglow flux densities to the burst physical parameters: 
the isotropic-equivalent kinetic energy (ZiK.iso), circumburst 
density (n), fractions of post-shock energy in radiating elec¬ 
trons (Ce) and magnetic fields (cb), and the electron power-law 
distribution index (p), with Nij) oc 7“^ for 7 > 7niin, where 
7min is the minimum Lorentz factor of the electron distribu¬ 
tion. 

The synchrotron spectrum is characterized by a flux nor¬ 
malization and three break frequencies: the self-absorption 
frequency (i/a), the peak frequency {v^), and the cooling fre¬ 
quency {vc)- Constraints on the physical parameters require 
knowledge of where the synchrotron break frequencies are 
located with respect to the observing bands. In most cases, 
there is not enough information to constrain the locations of 
Va and v,„ with respect to the observing bands, so we make 
assumptions about their locations (detailed in the next sec¬ 
tions). However, there are several cases in which we can use 
the available data to determine the location of Vc with respect 
to the X-ray and optical bands. To determine the location of 
Vc, we first determine the temporal and spectral power-law 
indices (a and /3, respectively, where iv oc from the 

X-ray and optical light curves and spectra. We then com¬ 
pare these indices to the standard relations given by the syn¬ 
chrotron model to determine whether Vc is located above or 
below the X-ray band. This also allows us to calculate the 
value of p, and governs how the flux es map to the burst phys¬ 
ical properties (iGranot & Sarill2003) . 

3.1. Temporal and Spectral Behavior 
3.1.1. X-rays 

To investigate the temporal behavior of the X-ray after¬ 
glows, we utilize ^^-minimization to fit a single power law 
model to each light curve in the form cx with tempo¬ 
ral index ax as the single free parameter and the best-fit flux 
normalization Co given by 


Co = 


^moddq'XFx,,' 

E V Cwdel.i 


( 1 ) 


where Cmodei,; are the un-normalized model fluxes, Fxj and 
ax.i are the observed fluxes and uncertainties, respectively, 
and N is the number of data points. Since early-time X- 
ray afterglow light curves are often subject to steep decays, 
plate aus, or flares which may contaminate the a fterglow emis- 
sion ( Nousek et al.l I2006t iZhang et"^ I2006t iMargutti et al.l 
l201lLl2013h , we on y utilize X-ray data at 5t > 1000 s, when 
bursts have typically settled into the power-law afterglow 
phase. For the X-ray light curves, we initially include all of 
the available data at 5t > 1000 s in the fit. In a few cases, 
there are light curve features beyond St ta 1000 s which sig¬ 
nificantly affect the fit: flares (GRBs 050724A and 111 121 A), 
plateaus (GRB 051221A), or steepenings (GRBs 051221A 
and 111020A). For these bursts, we exclude the time intervals 
that contain such features in the fits. In 10 cases, there only 
exists a single detection and an upper limit at St > 1000 s. 


Table 2 

Short GRB Spectral and Temporal Power-law Indices 


GRB 

ax 

Px 

<^opt 

/^opt 

050509B 

-1.10±0.25 

-0.88 ±0.34 



050709 

-1.23±0.10 

-1.24±0.35 

-1.42 + 0.08 


050724A 

-0.93 ± 0.08 

-0.81 ±0.15 

-1.74 + 0.11“ 

-0.82 + 0.03“ 

050813 

< -0.004 

-1 3+2-1 



051210 


- 2.1 ± 0.5 



051221A 

-1.08±0.12 

-1.0±0.2 

-0.97 + 0.06 


060121 

-1.23 ±0.20 

-1.07 ±0.16 

-0.60 + 0.24 


060313 

-1.47 ±0.39 

-0.96 ±0.09 

-0.70 + 0.19 

-1.35 + 0.19 

060502B 

060801 


-2 07+' '^*’ 

-0.68 ±0.12 



061006 

-0.85 ±0.30 

-0.90 ±0.25 

-0.43 + 0.08 


061201 

-1.96±1.18 

-0.66 ±0.12 



061210 

070429B 

-1.71±1.15 

<-1.12 

-1-86^:61 
-2.0 ±0.63 



070707 

<-0.65 


-2.55 + 0.22 


070714B 

-1.96 ±0.69 

-1.07 ±0.19 

-0.81+0.11 


070724A 

-0.95 ±0.33 

-0.60 ±0.25 

<-0.15 

-0.58 + 0.02 

070729 

070809 

071227 

-1.09±1.10 

-0.97 ± 0.27 

/T c+0.44 
^•'^-0.25 

-0 37+*’-*^ 

-0.90 ±0.31 

-0.73 + 0.33 

<-0.24 


080123 

-0.77 ±0.19 

-1.6±0.4 



080426 

-1.54 ±0.33 

-1.03±0.I6 



080503 

080702A 

<-0.38 

1 co+0.26 

i.jj_0.i2 

- 1.00 ± 0.43 



080905A 


- 0.53 ± 0 . 1 8 

-0.39 + 1.36 


080919 

081024A 

-1.23 ±0.69 

_i g±0.63 
^•^-0.38 

-0.70 zb 0.38 



081226A 

090305A 


9 97+1.24 
"^•^'-0.32 

-0.95 + 0.30 

-0.74 + 0.09 

-1.67 + 0.67'’ 
-0.71 + 0.27 

090426A 

-1.15±0.16 

-1.04 ±0.09 

-0.58 + 0.16 

-0.94 + 0.06* 

090510 

-1.79 ±0.63 

-0.75 ±0.08 

-2.37 + 0.29 

-0.85 + 0.05 

090515 

090607 

< -0.72 

1 co+0.78 
^•“^■^-0.28 
- 1.2 ± 0.4 

<-0.12 


090621B 

-1.48 ±0.54 

9 7+0.68 
^•'-1.0 



091109B 

-0.83 ±0.28 

- 1.1 ± 0.3 

-0.49 + 0.45 


100117A 


-1.6 + 0.3 

-1.60 + 0.33 


100206A 

100213 


-2 05+‘ -®^ 

^•^-^-0.53 
-2 0+^ -2 



100625A 

< 0.37 

-1.5±0.2 



100702A 


-1.7±0.3 



101219A 

-1.37±0.13 

-0.8 ±0.1 



101224A 

110112A 

-1.10±0.05 

—2 4‘*’1-^ 

-1.2±0.2 

<-0.32 


111020A 

-0.78 ±0.05 

-1.04 ±0.16 



111117A 

-1.21 ±0.05 

-1.0±0.2 



111121A 

-1.55 ±0.30 

-0.90 ±0.12 



111222A 

120305A 

120521A 

< -0.30 

-0 7+1+1 

^•'-0.4 

9 9+0.6 
^•^-0.4 

-0.8 + 0.25 



120630A 


-0.8 ±0.3 



120804A 

-1.02±0.10 

-1.1+0.1 



121226A 

-1.12±0.28 

-1.5 + 0.25 



130313A 

130515A 


-1 6+2* 

-0.7 + 0.31 



130603B 

-1.88±0.15 

-1.2 + 0.1 

-1.26 + 0.05 

-2.0 + 0.1* 

130716A 

< -0.52 

1 1+0.56 
-0.38 
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Table 2 — Continued 


GRB 

ax 

fix 

<^opt 

/^opt 

130822A 


-0.89^:‘4 



130912A 

-1.33±0.18 

-0.57 ±0.13 

<-0.42 


131004A 

-I.lib0.61 

-0.8 ±0.3 

-1.94±0.25 

-1.52±0.13* 

140129B 

-1.6±0.26 

-0.97 ±0.11 

-1.54^ 


140320A 

< -0.23 

o 0+2.6 
■^•°-0.88 



140516A 

-0.42 ±0.14 

-0.70 ±0.44 



140622A 

<-0.16 




140903A 

-1.05 ±0.20 

-0.60 ±0.10 



140930B 

-1.27±0.15 

-0.73 ±0.28 



141212A 


_1 7+2.75 
'-0.81 



150101A 

<-0.12 

-0.40 ±0.56 



150101B 

-1.07 ±0.15 

-0.67 ±0.17 

-1.01 ±0.62 


150120A 


-1.0±0.31 



150301A 


-0.7 ±0.13 




Note. — Error bars correspotrd to Ict confidence. When applicable, ax 
and Qopt represent pre-jet break values. Values of /3opt are observed values 
and are uncorrected for intrinsic rest-frame extinction. Ay. 

“ These values are computed over the same time interval as the X-ray flare 
that is super-imposed on the underlying afterglow power-law decay. No op¬ 
tical detections exist for the underlying afterglow. 

* This value corresponds to the spectral behavior after the jet break. 

No uncertainties are given for the optical light curve. 

SO we can only extract an upper limit for ax- The resulting 
best-fit values for ax, along with Icr uncertainties, are listed 
in Table|2] Also listed are the X-ray spectral indices, fix, from 
the relation fix = 1 -T. The X-ray afterglows have weighted 
mean values and Icr uncertainties of {ax) = -1.07;|;q3' and 
(/3x) = -0.89;!^'7g, and median values of ax,med ~ -1.16 and 

/5x,med ~ —1.06. 

3.1.2. Optical 

We determine the temporal index of the optical observations 
(ttopt, where/v.opt cx in the same manner as described in 
Section [3.1.1l using the hlter with the most well-sampled light 
curve for each burst. If there are multiple biters for which we 
can determine Oopt, we independently bt aopt for each biter 
and report the weighted mean (Table |2]l. A few short GRBs 
have well-measured optical spectral indices from contempo¬ 
raneous multi-band data, but do not have a well-sampled light 
curve in a single biter, preventing a measurement of ofopt. For 
these events, we use the measured value of /?opt (see below) to 
extrapolate all of the available afterglow data to a single bi¬ 
ter, and then determine the temporal decay index from these 
observations. In this manner, we are able to measure the opti¬ 
cal temporal decay index for 19 short GRBs, and place upper 
limits in bve cases for bursts with only a single detection and 
an upper limit. The best-bt optical temporal indices and Icr 
uncertainties are listed in Table |2l The weighted mean for the 
19 short GRBs with measured values is (oopt) = -1.07;|;Q g[, 


and the median is (Topt.med ~ -0.99. 

If there are contemporaneous observations in multiple bi¬ 
ters, we use these to determine the observed spectral slope, 
/3opt (Tv,opt oc For a few bursts, there are multi-band 

observations taken at different times, as well as a measure¬ 
ment of aopt from a well-sampled light curve in a single bi¬ 
ter; in such cases, we use aopt to interpolate the data from 
multiple biters to a common time. To determine /?opt, we 
then use x^-minimization to bt the optical/NIR photomeby 
to a power law model. The values of /3opt are listed in Ta¬ 
ble 12 The weighted mean for all short GRBs with deter¬ 
mined spectral indices after incorporating rest-frame extinc¬ 
tion (Section [3. 1.3b is (/3opt) = -0.71;|^-5j, and the median is 

/5opt,med ~ —0.88. 

3.1.3. Rest-Frame Extinction 

In four cases where there are contemporaneous observa¬ 
tions in multiple optical/NIR biters, a single power law model 
provides a poor bt to the broad-band photometry (x^ ^ 3). 
For these bursts, we include the line-of-sight rest-frame ex¬ 
tinction from the host galaxy (Ay°®'^) as a second free parame¬ 
ter in the bt. We c onstrain the extinctio n using the Milky Way 
extinction curve (ICardelli et al.lll98^ . and employ the burst 
redshifts (Table [T]i to obtain the rest-frame extinction. For 
bursts with no determined spectroscopic redshift, we assume 
z = 0. 5 set by the median of the short GRB population (iBerged 
l2014ft . We bnd non-zero extinction values for GRBs 060121, 
070724A, 081226A, and 130603B (Table |3]l. The resulting 
values for the optical spectral indices, uncertainties, and Ay”**' 
are listed in Tabled 

For bursts where we do not have enough information from 
the optical/NIR bands alone to constrain the specbal behavior 
of the afterglow, we initially assume that there is no rest-frame 
extinction. Ay”®* = 0. However, in eight cases, the difference in 
slope between the X-ray and optical bands is shallower than 
expected (e.g., |/3ox| ^ |/3x|-0.5, where fiox is the spectral 
index between the X-ray and optical bands), suggesting that 
either there is intrinsic extinction, or that the X-rays do not 
originate from the forward shock. Assuming the former ex¬ 
planation, we include the minimum amount of extinction re¬ 
quired until the X-ray and optical solutions agree to within the 
Icr uncertainties. In two cases, GRBs 080919 and 101219A, 
there is only an upper limit on the optical afterglow bright¬ 
ness; thus we can only determine a lower limit on the rest- 
frame extinction. We note that GRB 080919 has the highest 
value of rest-frame extinction, with Ay°'*'^ > 6 mag. However, 
this burst has a sightline close to the Galactic plane and there¬ 
fore has a highly uncertain Galactic extinction, which likely 
affects the inferred value for Ay°®l For the 12 events with 
rest-frame extinction, we bnd values of Ay°*' ss 0.3 -1.5 mag. 
These Ay°'*'^ values are also listed in Tabled 

3.2. Determination of Electron Power Law Index and the 
Location of the Cooling Frequency 
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Table 3 

Inferred Properties 


GRB 

/•host 


P 

es 

iso,52 

< £^K,iso,52 > 

Vj 

<n> 


(mag) 




(10^^ erg) 

(10^^ erg) 


(cm-3) 

050709 

0 

Y 

2.31 ±0.13 

0.1 

0.09 

2.6!5]j ^ 

0.97 

^•^-0.4 


0 

Y 

2.31 ±0.13 

0.01 

0.09 

6.2!0;^ X 10-3 

0.93 


050724A 

0 

N 

2.29±0.10 

10-^ 

0.24 

A 1 O+0.08 

^■^®-0.05 

0.58 

0.89!«:3« 

051221A 

0 

Y 

2.24 ±0.07 

0.1 

1.3 

^•^*^-0.01 

0.89 

A (-^->+0.006 

^•^”^-0.005 


0 

Y 

2.24 ±0.07 

0.01 

1.3 

^•^'-0.03 

0.83 

A ^4+0.05 
^•^^-0.04 

060121 

1.6 

Y 

2.24 ±0.20 

0.1 

4.5 

^•^^-0.04 

0.96 

5 . 4 + 1 * X 10-3 


1.6 

Y 

2.24 ±0.20 

0.01 

4.5 

^•^■^-0.04 

0.95 

^•^^-0.06 

060313“* 

0 

Y 

2.03 ±0.20 

0.1 

2.9 

A 4 C+O.O 5 

^•^^-0.05 

0.87 

3.3!3;0 X 10-3 

061006 

0 

N 

2.39±0.31 

0.1 

1.1 

A 54+0.85 
^•^^-0.37 

0.63 

2.2+}*, X 10-3 


0 

N 

2.39±0.31 

0.01 

1.1 

1 1+2-^ 

0.50 

1.2+fi X lO-* 

061201 

0 

N 

2.35 ±0.24 

0.1 

0.05 

A AC+0.10 
^•^^-0.03 

0.47 

5.0!**g X 10-3 


0 

N 

2.35 ±0.24 

0.01 

0.05 

^•^-0.09 

0.29 

2 . 7 ;™ X 10--* 

0707146*“ 

0.5 

Y 

2.30 ±0.35 

0.1 

1.7 

A 1+0.02 
'^•^-0.02 

0.94 

0 036+^-®^^ 
^•^^^-0.011 

070724A 

1.5 

N 

2.24 ±0.33 

0.1 

0.03 

'^•■^^-0.20 

0.07 

1 . 9 +pg X 10-3 


2.0 

N 

2.24 ±0.33 

0.01 

0.03 

1 1+3.0 

^•^-0.8 

0.02 

9.3pio X 10-3 

070809 

0 

N 

2.12±1.47‘' 

0.1 

0.09 

^•^-0.3 

0.14 

2.2+[3^ X 10-3 


0 

N 

2.12±1.47‘' 

0.01 

0.09 

1 1+2'? 

^•^-0.8 

0.07 

1.2+fi X lO-* 

071227 

0 

Y 

1.92±0.31 

0.1 

0.14 

8.4+2;5 X 10-3 

0.94 

1.9pp 


0 

Y 

1.92±0.31 

0.01 

0.14 

8.9!f;f X 10-3 

0.94 

60p5 

080426 

0 

Y 

2.29 ±0.26 

0.1 

0.35 

^•^‘^-0.01 

0.87 

0.04!P^ 


0 

Y 

2.29 ±0.26 

0.01 

0.35 

^•^'^-0.01 

0.85 

1 2+^-2 
^•^-0.5 

080905A 

0 

N 

2.06 ±0.36 

0.1 

0.02 

0.04^:‘3 

0.34 

i.3(;332 X 10--* 


0 

N 

2.06 ±0.36 

0.01 

0.02 

0.08^:;!^ 

0.21 

7.1!**“ X lO-* 

080919 

^6 

Y 

2.97 ±0.68 

0.1 

0.07 

^•^^^-0.004 

0.78 

0.20!“:*“ 



Y 

2.97 ±0.68 

0.01 

0.07 

u.uz,!^_0 006 

0.70 

5 - 11 :? 

081024A 

0 

N 

2.40 ±0.76 

0.1 

0.11 

A 1 1+0.22 
^•^^-0.07 

0.51 

8 . 1 ^ 37 “ X 10-3 


0 

N 

2.40 ±0.76 

0.01 

0.11 

^•^^-0.20 

0.31 

4.3!*“2 3< 10-* 

081226A 

1.0 

N 

2.27 ±0.39 

0.1 

0.09 

^•^^-0.13 

0.32 

3.2!f5 X 10-3 


1.0 

N 

2.27 ±0.39 

0.01 

0.09 

0 43+1-5 
^•^■^-0.33 

0.18 

1.7!34 X lO-* 

090426A 

0 

Y 

2.13±0.14 

0.1 

2.0 

1 4^-4 
^•^-0.3 

0.59 

A 04+0.04 
^•^^-0.02 


0 

Y 

2.13±0.14 

0.01 

2.0 

^•^-0.3 

0.57 

1 2+1-^ 

090510 

0 

N 

2.65 ±0.08 

0.1 

0.77 

A - 77 + 0.57 

^•''-0.33 

0.50 

i.2!3'5 X 10-3 


0 

N 

2.65 ±0.08 

0.01 

0.77 

1 9+3.0 
-1.2 

0.29 

6.4!P“ X 10-3 

090607 

0 

Y 

2.40 ±0.76 

0.1 

0.10 

2.1!“;^ X 10-3 

0.98 

0.84!«:32 


0 

Y 

2.40 ±0.76 

0.01 

0.10 

2.6!“;3 X 10-3 

0.98 


0906216“ 

2.5 

Y 

2.64 ±0.72 

0.1 

0.07 

U.UZJ>_0 006 

0.75 

^•^^-0.02 


2.5 

Y 

2.64 ±0.72 

0.01 

0.07 

0 03 

w.u.:)i_0 007 

0.68 

1 0+0-52 
-0.27 
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Table 3 — Continued 


GRB 

^host 

Uc < 

P 


iso,52 

< £^K,iso,52 > 

+ 

<n> 


(mag) 




(1033 erg) 

(1033 erg) 


(cm“^) 

091109B 

0 

N 

2.40 ±0.32 

0.1 

0.18 

o-25U:l^ 

0.42 

2.8+34 X 10-3 


0 

N 

2.40 ±0.32 

0.01 

0.18 


0.18 

9.4™ X 10-3 

100117A 

0 

Y 

2.36 ±0.30 

0.1 

0.22 

^•^^^-0.003 

0.92 

0 

^•^^-0.01 


0 

Y 

2.36 ±0.30 

0.01 

0.22 

u.uz:)_0 004 

0.90 

1 9+0.9 
^•^-0.3 

101219A 

^2.5 

N 

2.73±0.13 

0.1 

0.74 

0.3ti 

0.68 

4.6:t39 X 10-5 


^2.5 

N 

2.73±0.13 

0.01 

0.74 


0.46 

2.4+f3 X 10-4 

110112A'’ 

0 

Y 

2.49 ±0.07 

0.1 

0.03 


0.31 

X 10-3 

111020 A‘''’^ 

0.5 

Y 

2.08 ±0.32 

0.1 

0.17 

A ^0+0.09 

0.26 

4.5!:f;§ X 10-3 

111117A 

0 

Y 

2.27 ±0.07 

0.1 

0.55 


0.90 

8.3+^-^ X 10-3 


0 

Y 

2.27 ±0.07 

0.01 

0.55 


0.89 

A 25+0.16 

^•^2>_o.o7 

111121A 

N/A 

N 

2.87 ±0.21 

0.1 

2.1 


0.40 

8.2)|3^ X 10-s 


N/A 

N 

2.87 ±0.21 

0.01 

2.1 

Q Q+l 1.1 

0.20 

4.2^48^ X 10-3 

120804A^ 

2.5 

Y 

2.08 ±0.11 

0.1 

3.4 

1 i-fn.3 

^•^-0.2 

0.76 

3.2+34 JO-3 


2.5 

Y 

2.08 ±0.11 

0.01 

3.4 

rj ^+0.2 

0.60 

0 014^^'^^^ 
^•^^^-0.001 

121226A“* 

1 

Y 

2.50 ±0.37 

0.1 

0.37 

A .^±0.08 

0.37 

4.0!i;g 10^^ 

130603B 

1.2 

Y 

2.70 ±0.06 

0.1 

0.37 

A 1 1+0.02 
^•^^-0.01 

0.77 

A aa+0.04 
^•^^-0.03 


0.3 

Y 

2.70 ±0.06 

0.01 

0.37 

A ^5+0.02 
^•^^-0.02 

0.72 

A 3 ^+0.08 

^•■2-0.04 

130912A^ 

1.3 

N 

2.49 ±0.17 

0.1 

0.16 

u.:>u_o2o 

0.25 

1.7+” X 10-3 


1.3 

N 

2.49 ±0.17 

0.01 

0.16 

1 4+0-2 

0.11 

X 10-3 

131004A 

0 

N 

2.57 ±0.48 

0.1 

0.45 

1 9+0.56 
^•^-0.36 

0.28 

1.2+34 X 10-3 


0 

N 

2.57 ±0.48 

0.01 

0.45 

2-8^l5 

0.45 

6.5!‘‘2 X 10^' 

140129B 

0 

N 

3.00±0.19 

0.1 

0.07 

0.98^:*^ 

0.06 

3.Hf4 X 10-5 


0 

N 

3.00±0.19 

0.01 

0.07 

T q+12.0 

0.02 

1.6!^3^ X 10-4 

140516A 

0 

N 

2.40 ±0.88 

0.1 

0.02 


0.54 

1.0^34^5 X 10-4 


0 

N 

2.40 ±0.88 

0.01 

0.02 

A a4_+0.18 

^•^^-0.03 

0.34 

5.5!^‘4 X 10-4 

140622A 

0 

N 

2.10±0.60 

0.1 

0.07 

A 1 2+0-26 
^•^^-0.08 

0.36 

5.8(;|* X 10-3 


0 

N 

2.10±0.60 

0.01 

0.07 

A 2S+0-97 

^•^2>_o20 

0.21 

3.2+‘'+3 X 10-4 

140903A 

0 

N 

2.27±0.16 

10-3 

0.08 

2-9t?I 

0.03 

3.4+34 X 10-3 

140930B 

0 

N 

2.67±0.19 

0.1 

0.40 

0.28^:« 

0.58 

4.6!^*3 X 10-3 


0 

N 

2.67±0.19 

0.01 

0.40 

1 8+1-7 
^•°-0.9 

0.18 

1.8iJ5 X 10-3 

150101B'' 

0.5 

N 

2.40 ±0.17 

0.1 

4.0 X 10*3 


6.6 xlO-3 

8.0(|4j X 10-s 


0.5 

N 

2.40 ±0.17 

0.1 

4.0 X10-3 

3.0!!:] 

1.3 xlO-3 

5.3+‘'5 X 10-6 


Note. — Quoted uncertainties are Icr. All solutions presented here are 
for a fixed = 0.1 and assume a lower density bound of = 10 “^ cm“^. 

^ We assume a redshift of ^ = 1 for this burst. 

^ No valid solution is found for eg = 0.01. 

^ Value of is determined from a comparison of the optical and X-ray 
bands, and not directly from the optical/NIR SED. 

^ Determined from ax alone. 


To determine the electron power law index /?, we use ax 
and px to constrain the location of the cooling frequency, 
with respect to the X -ray band. We use the relations given by 
iGranot & Saril (I2QQ2I) which relate ax and Px to the value of 
p for the scenarios <'t^x <'^c and Uc < t^x by 
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Figure 2. Histogram of electron power law index, p, for 38 short GRBs as 
inferred from the X-ray temporal and spectral indices. The weighted mean 
and la uncertainties for the population (blue arrow) is {p) = 2A3'^^. These 
values correspond to those listed in TableO 


r l-2/3z] 

p = { A V v^<vx<i'c (2) 

f -2&1 

p = < l-^Xax > Vc< vx- (3) 

For a given burst, we calculate the values of p and Icr uncer¬ 
tainties using Equations |2] and [3 and standard propagation of 
errors. We select the valid scenario under the condition that 
the values of p independently determined from ax and jSx for 
a given scenario agree within the Icr uncertainties. Following 
this condition, we can constrain the location of Vc with respect 
to the X-ray band for 38 bursts (TableO. Using Equations |2] 
and[3l we calculate the weighted mean for the value of p in 
the valid scenario; the resulting values and uncertainties are 
listed in Tableland displayed in Eigure|2l We note that in the 
case of GRB 071227, the condition is satisfied for Vc < vx, but 
has a median of p = 1.92 ± 0.31, which yields a divergent total 
integrated energy. Thus for this burst, we employ p = 2.05 in 
our subsequent analysis (Table[3ll. 

Under the reasonable assumption that the optical band lies 
between the peak frequency, v^, and the cooling frequency, 
Vc (i.e., Vm < z/opt < Vc), we use the available values for Oopt 
and /3opt and Equation|2]to independently determine the value 
of p. We then include this in our weighted average of p for 
each burst (Eigure|2]i. The weighted mean for the sample of 
38 bursts is (p) = 2.43;|;o28 Eigure|2]i. 

Eor the remaining bursts in the afterglow catalog, there is 
not enough information to determine the location of the cool¬ 
ing frequency with respect to the X-ray band, and therefore 
the correct value of p. We thus concentrate on the subset of 
38 bursts with determined values of p for our subsequent anal¬ 
ysis. We find Vc < vx for 18 cases, while Vc > vx for 20 cases. 

3.3. The Isotropic-Equivalent Kinetic Energies and Densities 
for Individual Bursts 


In t he standard synchrotron model from iGranot & S aril 
(l2002h . the dependencies on the isotropic-equivalent kinetic 
energy, circumburst density, and the microphysical parame¬ 
ters are as follows for a given flux density, F^ i and observing 
band, vg. 


oc 


1/2 5/6 -2/3 1/3 

"O ^K.iso,52'^c-l '^B-1 

n^'^p 4 p-l 4 

^K,iso,52^^,-1 ^5,-1 

P 4 fP ^ f i 

'^K,iso,52'^e,-UB,-l 


Va < Vi < Vm 
Vm < Vi < Vc 
Vi > Vc 


(4) 


where no is in units of cm“^, £’K.iso.52 is in units of 10^^ erg, 
and Ce-i and es.-i are in units of 0.1. In addition to these four 
parameters, F^^i is also dependent on the redshift, luminosity 
distance, 12,, time after the burst, 5 t, and the value of p; th e ex¬ 
act dependencies are provided in iGranot & Saril (120021) . We 
note that for 12, > Vc, the flux density is independent of cir¬ 
cumburst density. Eor bursts with no spectroscopic redshift, 
we assume z = 0.5. In all cases, we cannot independently con¬ 
strain Cg and cb since this requires knowledge of the locations 
of the three break frequencies {va, v^, and Vc), which gen¬ 
erally necessitates well-sampled light curves and spectra in 
multiple bands. Thus, in order to determine ranges for Ek.Iso 
and «, we fix the values of the microphysical parameters. We 
first consider the fiducial case that = 0.1 and = 0.1. 

Eor each burst, we determine the constraints on Ek/so and 
n by computing individual probability distributions for each 
observation. We then assign the probabilities to a grid of 
values, and use a joint probability analysis to calculate the 
distributions in each parameter. Eor the grid, the ranges of 
the density and isotropic-equivalent kinetic energy are n, = 
10“®-10^ cm“^ and £,■ = 10"^®-10®"^ erg, with 1000 logarithmi¬ 
cally, uniformly-spaced steps in each parameter. We choose 
the lower bound of the density range, = 10“® cm“^ to 
match the typical density of the intergalactic medium (IGM). 


3.3.1. Detections 

To calculate the individual probability distributions for af¬ 
terglow detections, we apply Equation |4] to the observations 
using the relevant regime for each observing band, Z2,. Eor the 
X-r ay ban d, we use the location of Vc as determined in Sec¬ 
tion [3TT] to determine which branch of Equation |4] to use. 
Since the value of p is primarily determined from the X-ray 
band, the Ek.Iso “ n relation remains unchanged when using 
different X-ray observations that follow the same temporal 
decline; thus, we only use one X-ray observation per burst. 
Eor the optical band, we make the reasonable assumption that 
Vm < Vopt < Vc and utilize the second branch of Equation |4l 
In some cases, individual optical/NIR observations lead to 
Ex.iso “ n relations that do not overlap within their 1 a uncer¬ 
tainties. In these cases, we use the weighted mean and stan¬ 
dard deviation of these relations (i.e., systematic uncertainty) 
as the opticaENIR solution. Eor the radio band, we assume 
that Va < z2iadio < v„t, (flrst branch of Equation [Hi to calculate 
the £K,iso“« relations. Since the radio band lies on a different 
segment of the synchrotron spectrum than the other bands, it 
contributes a Ek/so-w relation with a different slope, and thus 
enables tighter constraints on the physical parameters (Eig- 
ureO. 

After calculating the unique probability distribution from 
each of the radio, opticaENIR and X-ray bands, we normalize 
the area under each of the distributions to unity. We assume 
that the uncertainties in the flux densities are Gaussian, and 
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Figure 3. Isotropic-equivalent kinetic energy versus circumburst density for four short GRBs with radio afterglow detections. In each panel, the X-rays (light 
blue), optical (orange) and radio (red) provide independent constraints on the parameter space. Measurements are shown as solid regions, where the width of 
the region corresponds to the Icr uncertainty. Upper (lower) limits set by the cooling frequency are denoted by blue dotted (dash-dotted) lines. The regions of 
parameter space ruled out by the observations are denoted (grey hatched regions). The median solution and Itr uncertainty is indicated by the black cross in each 
panel, corresponding to the values listed in Table[^ For each burst, the joint probability distributions in n (bottom panel) and iiK.iso (right panel) are shown. Red 
lines correspond to the median, and dotted lines are the la uncertainty about the median. The green line corresponds to £-y,iso. The only optical observations 
available for GRB 050724A are during a flare; thus we do not include the optical data in our analysis. The addition of a detection in the radio band is crucial in 
constraining the best-fit solution, and in three cases constrains eg. 


thus each band contributes a unique, log-normal distribution. 
These distributions are shown for the four bursts with radio 
afterglow detections (Figure [3), and 34 bursts with radio af¬ 
terglow non-detections: 16 bursts with i/c < vx (Figure|4|i and 
18 bursts with Vc > vx (Figure|5]l. 

3.3.2. Upper Limits 

For flux density upper limits, we use Equation |4] to deter¬ 
mine the Ej^ iso-n relation at the 3cr upper limit. We then 
assign zero probability to the £’K,iso“n parameter space above 
the relationship, and assign a constant probability to the al¬ 
lowed parameter space below the relationship, normalized to 
unity. The upper limits are denoted by dashed lines in Fig¬ 
ures [311] where the colors of the lines correspond to the rel¬ 
evant observing band (X-rays: blue, optical: orange, radio: 
red). Regions of parameter space that have been ruled out by 
the upper limits are marked as hatched regions. 

3.3.3. Cooling Frequency Constraint 

We utilize the relative location of the cooling frequency as 
a hnal constraint, since it depends on a combination of energy 


and density, by 

Vc OC tlQ £K,iso,52^B,4 > (5) 

with addition al _ dependencies on 6t and redshift 

(iGranot & S aril |2003)- For the cases in which the X-ray 
band is located above the cooling frequency (vc < vx), we 
employ a maximum value at the lower edge of the X-ray 
band, r'e,max = 7.3 x 10'® Hz (0.3 keV), to obtain a lower limit 
on the combination of energy and density. The corresponding 
probability distribution has zero value in the param¬ 

eter space below the relation, and a constant value above the 
relation, where the area in the allowed parameter space is 
normalized to unity. The lower limits are shown as as blue 
dot-dashed lines for GRBs 051221A and 130603B in Figure|3 
and for all bursts in Figure |4| In cases where the X-ray band 
is below the cooling frequency {v^ < vx < Vc), we set the 
cooling break to a minimum value, i/c,mm = 2.4 x 10*^ Hz 
(10 keV) at the upper edge of the X-ray band, and determine 
the £’K,iso “ n relation for each burst using Equation [5] This 
constraint sets an upper limit on the combination of energy 
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Figure 4. Isotropic-equivalent kinetic energy versus circumburst density for 16 short GRBs with solutions for Vc < I'x assuming fiducial values for the micro¬ 
physical parameters of = ea = 0.1. In each panel, the X-rays (light blue), optical (orange) and radio (red) provide independent consteaints on the parameter 
space. In particular, the X-ray band provides an estimate of iiK.iso- Measurements are shown as solid regions, where the width of the region comesponds to the la 
uncertainty, while upper limits are denoted as dashed lines. Setting the cooling frequency to a maximum value of t'c,max = V.2 X 10*® Hz (0.3 keV) provides an 
additional constraint (dark blue dot-dashed line). The regions of parameter space ruled out by the observations are denoted (grey hatched regions). The median 
solution and Icr uncertainty is indicated by the black cross in each panel, corresponding to the values listed in Table For each burst, the joint probability 
distributions in n (bottom panel) and Fk.Iso (right panel) are shown. Red lines correspond to the median, and dotted lines are the 1 a uncertainty about the median. 
The green line corresponds to £-,,iso. 


and density. We form the probability distribution in the same 
manner as for afterglow upper limits. The limits for each 
burst set by the cooling frequency are shown as blue dashed 
lines for GRBs 050724A and 140903A in Figure [3 and all 
bursts in Figure |5] 

3.4. Joint Probability Distributions 

Since each of the observing bands, as well as the location of 
the cooling frequency, contribute an independent probability 
distribution, we calculate the joint probability from a product 
of these distributions for each burst. To obtain 1-dimensional 
probability distributions in isK.iso and n, we integrate over each 
of the parameters. Finally, we normalize the area under each 
1-dimensional distribution to unity. The resulting distribu¬ 
tions, Pin) and FCfiK.iso). for 34 bursts are shown in Figures|4]- 
|5]for the fiducial microphysical parameters, = 0.1. The 

median values and Icr uncertainties in isotropic-equivalent ki¬ 
netic energy and circumburst density are also shown in these 
figures and listed in Table[3] 

Figure [2 shows the probability distributions for the four 
events with radio afterglow detections. In three of the four 


cases, we can use the available data to place additional con¬ 
straints on tB, fixing = 0.1. For GRBs050724A and 
140903A, the afterglow data require that eg < 10“"^ and < 
10“^, respectively. For larger values of e^, the constraint from 
the cooling frequency becomes more stringent and conflicts 
with the solutions obtained from the afterglow observations 
on the grid of allowed values. For GRB 051221 A, we find a 
significantly better fit at = 0.01. In all cases, the addition of 
the radio band enables tighter constraints on both the energy 
and the density. 

4. DENSITY AND ENERGY SCALE FOR SHORT GRBS 

To quantify the distributions of circumburst densities and 
isotropic-equivalent kinetic energies for the entire sample, 
we calculate the combined probability distributions from the 
sum of the 1-dimensional probability distributions. Pin) and 
PjiiK.iso)- We sum the individual distributions for all bursts 
with valid solutions at eg = 0.1, as well as GRB050724A 
icB = 10“^) and GRB 140903A (eg = 10“^), to create cumula¬ 
tive probability distributions for both density and kinetic en¬ 
ergy, shown in Figure |6] 
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Figure 5. Isotropic-equivalent kinetic energy versus circumburst density for 18 short GRBs with solutions for Vc > assuming fiducial values for the micro¬ 
physical parameters of Cc = eg = 0.1. In each panel, the X-rays (light blue), optical (orange) and radio (red) provide independent consteaints on the parameter 
space. Measurements are shown as solid regions, where the width of the region corresponds to the Icr uncertainty, while upper limits are denoted as dashed lines. 
Setting the cooling frequency to a minimum value of 1^^,111111 = 2.4 X 10** Hz (10 keV) provides an additional constraint (dark blue dashed line). The regions of 
parameter space ruled out hy the observations are denoted (grey hatched regions). The median solution and Itr uncertainty is indicated by the black cross in each 
panel, corresponding to the values listed in Table[3] For each burst, the joint probability distributions in n, with an imposed lower bound of «ini„ = 10^ cm“* 
(bottom panel), and isi, (right panel) are shown. Red lines correspond to the median, and dotted lines are the la uncertainty about the median. The green line 
corresponds to £-,,,150 ■ 


The median values for the population of 38 bursts are 
(n) « 2.9 X 10“^ cm“^ and (iiK.iso) « 1.9 x 10^^ erg (Figure| 6 ] 
and Table m. The density which corresponds to 90% of the 
cumulative distribution (ngo) is « 0.4 cm“^, and the fraction 
of probability that lies at densities of < 1 cm“^ (fn<icm-^) is 
0.95 (Figure| 6 ]and Table|4]l. We also calculate these statistics 
for the sub-sample of bursts with < ixx, corresponding to 
the events in Figure 0] that utilize all three branches of Equa¬ 
tion 0 ] and therefore have relatively well-constrained energies 
and densities (compared to the broader pro babilit y distribu¬ 
tions for events with zz^ > zz^; see Section 0.2.3b . We find 
that this distribution has a median of (zz) r;;! 3.7 x 10“^ cm“^. 


zzgo « 1.0 cm and that /„<icm -3 ~ 0.90. Furthermore, the 
median value for kinetic energy decreases by a factor of 2 , 
to Ek.Iso ~ 9.6 X 10^° erg. The cumulative distributions in 
density and energy are shown in Figure | 6 ] and the population 
statistics are listed in Table 0] 

4.1. Isotropic-equivalent ■j-ray energy and Efficiency 

We compare the inferred isotropic-equivalent kinetic en¬ 
ergy to the 7 -ray energy by computing the isotropic- 
equivalent 7 -ray energy, to represent the energy range 

Ri 1-10"^ keV (to match the widest energy ranges for current 
GRB detection satellites). 
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Figure 6. Combined and cumulative probability distributions in n and £’K,iso assuming Se = es = 0.1 (top) and €e = 0.1, es = 0.01 (bottom). All scenarios also 
include GRBs050724A and 140903A, with eb = 10“^ and 10“^, respectively. Each panel shows three populations: all bursts with an imposed lower bound of 
njnin = 10“^ cm“^ (red), all bursts with an imposed lower bound of «min = 10^ cm“^ (black dot-dashed), and the sub-sample of events with Uc < I'x (light grey 
dashed). Color-coded an'ows from the bottom denote the median for each distribution (and in some cases, staggered for the sake of clarity), and lines denote 90% 
upper limits. For = 0.1 (e^ = 0.01), the cumulative distributions indicate that for w > 3 x 10“^ cm“^ {n > 0.3 cm“^), the distributions are virtually independent 
of the choice of provided that > 10^ cm“^. This allows us to place robust 90% upper limits of n < 0.5- 1.0 cm“^ (« < 3-23 cm“^) for our sample 
of bursts. In addition, ^ 95% (~ 80%) of the total probability for all events lies below densities of 1 cm“^. The range of median isotropic-equivalent kinetic 
energies for all scenaiios is ^ (2-20) x 10^^ erg (c.f., TableHJ. 


Andl 

» “ ^bol ^ 

1 +z 


erg 


(6) 


where A:boi is the bolometric correction factor to convert the 
fluence to an energy range of ~ 1 - lO'^ keV, di is the lu¬ 
minosity distance in cm, and is the fluence in units of 
erg cm“^. For cases in which the fluence is calculated over 
the 15-150 keV Swift energy range, we use kboi = 5. If a 
burst is detected by other 7 -ray satellites which cover a wider 
energy range of ss 10- 1000 keV (e.g., Fermi, Konus-Wind, 
Suzaku), we utilize the measured fluences from these satel¬ 
lites and kboi = 1 to calculate the 7 -ray energy. Individual 
values of FC-y iso are shown in Figures [3^ and listed in Ta¬ 
ble |3l while the distribution for 38 bursts is shown in Fig¬ 
ure IT). We find that the range in isotropic-equivalent 7 -ray 
energy is (0.04-45) x 10^' erg, with a median value of 
{E^^iso) ~ 1.8 X 10^* erg, similar to the ranges and median 


values of the isotropic-equivalent kinetic energy (Figure|2ll 
We also calculate the 7 -ray efficiency. 


Vj = 


^7,ISO 


-'7,ISO 


+ EK.i 


(7) 


as well as the la uncertainties in rjj for each burst, follow¬ 
ing standard propagation of errors from the Icr uncertainties 
in ZiK.iso- The resulting values of are listed in Tableland 
the cumulative distributions, after incorporating the Icr uncer¬ 
tainties, are shown in Figure [T] We find a wide range in r],y, 
Ri 10“^- 1, and note that the lower bound is set by the sin¬ 
gle outlier, GRB 150101B (TableO. Excluding this burst, the 
lower bound is « 0.03 (GRB 140622A). 

The distribution for = 0.1 has a median of =0.56^g 3 ® 
(Icr uncertainties). The isotropic-equivalent kinetic energy 
scale for the eg = 0.01 case is comparatively high (c.f.. Fig¬ 
ure 01 , and thus the median value for the 7 -ray efficiency is 





















Short GRB Afterglows 


15 



Figure 7. Left: Histograms of isotropic-equivalent 7-ray energy (red), and infen'ed isotropic-equivalent kinetic energies assuming es = 0.1 (dark grey, thick bars) 
and eg = 0.01 (light grey, thin bai's) for 38 short GRBs from broad-band afterglow modeling. Median values for each population are denoted by the color-coded 
arrows from the top. Values for the isotropic-equivalent energies correspond to those listed in Table[3] Right: Cumulative distributions of 7-ray efficiency 
assuming = 0.1 (dark grey) and eg = 0.01 (light grey), taking into account the \a uncertainties in for each burst. Median values are listed in the figure and 
denoted by the color-coded arrows from the bottom. Arrows are staggered for clarity. 


relatively low, ( 77 ,,,) = 0. 40^035 (Icr uncertainties; Figure]?]). 

4.2. Alternative Cases 

We have thus far made assumptions about the values of the 
microphysical parameters (e^ = = 0.1), the redshifts (z = 0.5 

unless otherwise determined), the minimum allowed density 
(Hmin = 10“® cm“^), and the cooling frequency (0.3 keV for 
t'c < vx', 10 keV for Vc > vx)- To explore the impact of these 
assumptions on our resulting distributions for the kinetic en¬ 
ergy and circumburst density, we consider alternative values 
for these parameters. 

4.2.1. The Value of es 

In some cases, a valid solution using the fiducial input of 
cb = 0.1 (with fixed Ce = 0-1) cannot be found. For instance, 
three of the four bursts with radio afterglows require that 
cb < 0.1 (FigurejSjand TablejSj). At a fixed value of = 0.1, 
cb is constrained to < 10“'^ for GRB050724A, < 10“^ for 
GRB 051221 A, and < 5 x lO'^ for GRB 140903A (Figure|3]). 
For the remaining bursts, we do not have enough informa¬ 
tion to rule out the = 0.1 scenario. Thus, we consider two 
additional cases for all bursts; cb = 0.01 and = 10 ““* (with 
fixed Ce = 0.1). Of the 38 bursts in our sample, 33 have valid 
solutions for = 0 . 01 ; the median and uncertainties in cir¬ 
cumburst density and kinetic energy for each burst are listed 
in Table |3| To create cumulative probability distributions, we 
repeat the same exercise as described in Section|4|for the 33 
events with valid solutions, displayed in Figure | 6 | 

Using eB = lO"”* as the fiducial value, the constraint from the 
cooling frequency conflicts with the solutions allowed by the 
afterglow observations in 16 cases, indicating that data do not 
allow such a low value of cb for these bursts. The population 
statistics for the 22 bursts with valid solutions are listed in 
Table H 


population with known redshifts (|Berge3|20T3)- However, in 
three cases, GRBs 060313, 111020A, and 121226A, we do 
not And a valid joint solution at z = 0.5. We And that there are 
valid solutions at redshifts of z > 1 , and thus assume z = 1 for 
these bursts. 

4.2.3. The Value o/nmin 

In the 20 cases in which the cooling frequency lies above 
the X-ray band { 1 /^ < vx < Vc), the X-ray and opticaFNIR 
bands occupy the same spectral regime (branch 2 of Equa¬ 
tion S and the resulting relations have the same 

slope (Figure [3] and |5]). Thus, the lower bound on the den¬ 
sity is set by our minimum grid value of Wmin = 10 “^ cm“^; 
the density is otherwise unconstrained at the low end and re¬ 
sults in broad probability distributions in both circumburst 
density and energy. To understand the impact of our choice 
of Hniin on the resulting distributions, we repeat the individ¬ 
ual probability analysis, employing a more stringent lower 
bound of Hmiri = 10 “"^ cm“^, at the low end of gas densities 
for the interstellp me dium (ISM; iKomi et ^ 119991; iMuralil 
1200(1 iGent et al.ll2013h . Since kinetic energy and density are 
inversely related, the upper bound on Ex.iso for each burst is 
naturally set by our choice of n^nm- We consider this alterna¬ 
tive value of Mmin for = 0.1 and = 0 . 01 . 

4.2.4. The Value of Vc 

In all cases, we have assumed that the cooling frequency is 
on the edge of the X-ray band (0.3 keV for Vc < vx', 10 keV 
for Vc > Vx). To test whether this assumption has any im¬ 
pact on our results, we repeat the individual and joint prob¬ 
ability analysis assuming that the cooling frequency is at 
the logarithmic mean of the 0.3- 10 keV Swift X-ray band, 
Vc.mid = 4.1 X 10 '^ Hz (1.7 keV). We consider this alternative 
value of the cooling frequency for = 0.1 and = 0 . 01 . 


4.2.2. Redshift 


4.2.5. Trends 


For bursts with no determined spectroscopic redshift, we 
have assumed z = 0.5, set by the median of the short GRB 


Taking these alternative values into account, we repeat the 
individual and joint probability analysis for each burst for nine 
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Table 4 

Circumburst Density and Kinetic Energy Population Statistics 


Scenario 

No. of Events 

{n) 

W90 

fb 

J n<l cm“3 

(£K,iso) 



(cm-3) 

(cm-') 


(erg) 

eB = 0.1 

All bursts, Wniin = 10“® 

38 

2.9 X 10^^ 

0.40 

0.95 

1.9 X 10'* 

All bursts, Wmin = lO^"* 

37 

3.7 X 10^^ 

0.49 

0.95 

1.2 X 10'* 

All bursts, tiuiin = 10“® cm“^, Vc= 1.7 keV 

38 

2.2 X 10^^ 

0.20 

0.95 

1.6 X 10'* 

Bursts with Vc < vx 

19 

3.7 X 10^^ 

1.01 

0.90 

9.6 X 10“ 

eB=0.01 

All bursts, Wmin = 10^® 

33 

5.2 X 10^^ 

2.7 

0.79 

2.9 X 10'* 

All bursts, ticXm = 10^ 

32 

1.4 X 10^^ 

3.1 

0.78 

1.7 X 10'* 

All bursts, Wmin = 10“® cm'^, Vc=l.l keV 

34 

1.6 X 10^^ 

2.1 

0.83 

2.4 X 10'* 

Bursts with Vc < vx 

14 

0.96 

23.3 

0.51 

6.4 X 10'“ 

es = 1 X 10^'* 

All bursts, Wmin = 10 ® cm ^ 

22 

3.0 X 10^^ 

771 

0.68 

1.8 X 10'' 


Note. — All scenarios include GRBs050724A and 140903A, 
with efl = 10^ and 10“^, respectively. 

“ This is the circumburst density which corresponds to 90% of the 
cumulative distribution. 

This is the fraction of the circumburst density cumulative distribu¬ 
tion below a value of 1 cm“^. 

different sets of input parameters in total. The population me¬ 
dians, as well as values for ngo and /,j<icm-3 are listed in Ta¬ 
ble |4l In addition, cumulative distributions for kinetic energy 
and circumburst density for varying values of eg and Wmin are 
shown in Figure |6] 

Overall, we find that a change in eg results in an increase 
in the circumburst density. For instance, assuming Wmin = 
10 “^ cm“^ for all bursts, decreasing by a factor of 10 to 
cb = 0.01 results in an increase in the median density by a 
factor of ss 1.8. This trend becomes more drastic for other 
scenarios; when assuming Wmin = 10“^^ cm“^, the median den¬ 
sity increases by a factor of « 3.8, and when considering only 
the bursts with Vc < vx, the factor of increase is « 26. When 
comparing the values of ngo and /„<icm-3, we also find overall 
trends commensurate with an increase in circumburst density. 
In particular, the factor of ten decrease in causes /„<icm-3 
to decrease by a larger factor, from « 0.95 to Ri 0.8 for all 
bursts. When considering only bursts with Vc <vx, fn<\cm-^ 
decreases from « 0.90 to Ri 0.51. The effect of on the me¬ 
dian kinetic energy is less pronounced, with increases by fac¬ 
tors of R:! 1.5 in all cases except when considering only bursts 
with Vc <vx', in this case, the median undergoes a slight de¬ 
crease by a factor of Ri 1.5 (Table |4|i. Overall, we find that a 
decrease from = 0.1 to 0.01 results in higher median den¬ 
sities by factors of Ri 2-25 depending on the considered sce¬ 
nario, and a uniform increase in the density cumulative distri¬ 
butions (Figure |6]l. When considering the more extreme input 
of cb = 10“^ for all bursts, the median density and kinetic en¬ 
ergy both increase by factors of R:^ 10, compared to the = 0.1 
case (Table IDl. 

We next investigate the effects of the values of Wmin and Vc 
on the parameter distributions. We find that at constant e^, 
the upper r:! 30-50% of the density cumulative distributions 
are virtually independent of our choice of Wmin, provided that 


fimin ^ 10 ^ cm ^ (Figure |6ll. Importantly, this also allows 
us to place robust 90% upper limits on the density that are 
largely unaffected by our choice of rijam- ngo ~ 0.4-0.5 cm“^ 
for cb = 0.1 and ngo R:! 2.7-3.1 cm“^ for cb = 0.01 (TablelHi. 
We find that the choice of cooling frequency within the X-ray 
band has a minor effect on the circumburst density, either a 
R:^ 1.3 factor of decrease {es = 0.1) or a factor of Ra 3 increase 
{eB = 0.01), while the median kinetic energy only decreases 
by a factor of Ra 1.2 in both cases (Table|4|i. 

When considering all cases with less extreme values of cb 
for all bursts, the median density range is (3-15) x 10“^ cm“^, 
with 90% upper limits of ngo Ra 0.4-3 cm“^. Furthermore, 
Ra80-95% of the probability is below Ra 1 cm“^ regard¬ 
less of the input parameters considered (Table |4|i. The me¬ 
dian isotropic-equivalent kinetic energy ranges from Ri (1.2- 
2.9) X 10^' erg. Including more extreme scenarios like the 
subset of bursts with Vc < vx and es = 10“^, the full median 
density range is Ra (3- 1000) x 10“^ cm“^ and the median 
isotropic-equivalent kinetic energy range is Ra (0.6 - 20) x 
10^2 erg (TableH. 

5. DISCUSSION AND IMPLICATIONS 
5.1. Afterglow Properties 

The X-ray, optical, and radio afterglow catalogs (Ta¬ 
bles A1-A3) allow us to analyze the observational after¬ 
glow properties of short GRBs as a population. By fitting 
the afterglow light curves, we find a weighted mean pre¬ 
jet break decline rate of {ax) ~ -1.07 at 5t > 1000 s for 
bursts with measured temporal indices, similar to the pre- 
j et break declines measu r ed from long GRB l ight curves 
(iNvsewander et all l2009t iRacusin et alJ 120091 iKann et ^ 
l2010t IZaninoni et ai.l 120131) , and slightly shallower than the 
value of ay Ki -1.2 found for 11 short GRBs in an earlier 
study (INvsewander et al.]l2009^ . We measure the optical de¬ 
cline rates and find the same weighted mean decline rate of 
(ciopt) ~ -1.07 from 19 well-sampled bursts with measured 
indices. 

From spectral fitting of the optical, near-IR, and X-ray 
data, we find 12 short GRBs which require rest-frame extinc¬ 
tion, with measured values of r:^ 0.3 -1.5 mag, and two 
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bursts with lower limits of >2.5-6 mag (TableO. We 
note that GRB 080919 has the highest value of rest-frame ex¬ 
tinction; however, the sightline to this burst is close to the 
Galactic plane and therefore has a highly uncertain Galac¬ 
tic extinction which likely affects the measurement of 
Prior to this study, evidence for > 0.5 mag of extinction 
has only been repor t ed in three shor t GRB s: GRB 070724A 
(iBerger et all l200^ iKocevski et all |2010h, GRB111020A 
jpong et al.112012^ and GRB 120804A dBerger et akllMl . 
Afterglow modeling in this work results in the same con¬ 
clusions for those three events, and includes nine additional 
events with rest-frame extinction. We note that most of 
the events with non-zero extinction and robust host associ¬ 
ations are in star-forming host galaxies; the single excep¬ 
tion is GRB 150101B which has an inferred value of Ay°*‘ 
0.5 mag and is located on the outskirts of an early-type galaxy 
(Fong et al., in prep). In comparison, « 15-20% of Swift long 
GRBs have optically sub-luminous afterglows that have been 
attributed to dust extinction. For long GRBs, inferred val¬ 
ues of Ay°‘*‘ ss 0.5 mag are common, with a substantial frac¬ 
tion of events wit h Ai*°^‘ > 1-2.5 mag (ICenko et alJ 120091 ; 
iPerlev et al.lf2()()^[2?illi . 

While rest-frame extinction can be explained by dust in 
star-forming regions in the local environments of long GRBs, 
substantial extinction in short GRBs cannot be easily ex¬ 
plained in the context of the double compact object merger 
progenitor. It is possible that these events are “impostors” 
which in actuality have massive star progenitors; in that case, 
we might expect them to be distinct in their 7 -ray proper¬ 
ties with longer durations and softer 7 -ray spectra. How¬ 
ever, there does not appear to be a correlation between short 
GRBs with extinction and their durations as they span the 
full range, with Tgo « 0 . 02-2 sec, and only three events 
with non-zero extinction have Tgo > 1 sec (GRBs 060121, 
070714B, and 121226 A). Furthermore, a study conducted by 
iBromberg et al.l (1201 3h assigned six of these events probabil¬ 
ities of having non-collapsar progenitors based on their 7 -ray 
properties. According to this study, 4/6 events have > 60% 
probabilities that they do not originate from collapsars. Thus, 
it is unlikely that the majority of these events are in fact “im¬ 
postors”. Instead, this suggests that some short GRBs may 
originate in a star-forming regions, or have progenitor sys¬ 
tems that can produce dust. 

5.2. Opening Angles 

Most well-sampled short GRBs exhibit a single after¬ 
glow decline rate in the X-ray (at 5t > 1000 sec) and 
optical bands. However, there are four short GRBs 
(GRBs051221A, 090426A, 111020A, and 130603B) which 
have temporal steepenings on timesc ales of 5t k, 2- 5 days , 
attributed to jet breaks (Tab l e fSl ISoderberg et al.l l2006l 
iNicuesa Guelbenzu et al.ll201 U iFcmg et al.ll2012^bl |2014|) . Jet 

breaks are achromatic features, and in principle can be de¬ 
tected in the X-ray through radio bands. The measurement 
of jet break time, in conjunction with the energy, density, and 
redshift, ca n be converted to a j et opening angle, 0 j, using 
(iSari et al.lll999HFrail et a]M)m 

9j = 9.5 t%\l deg, ( 8 ) 

where tjo is in days. The opening angle measurements for 
the four short GRBs with jet break detections are listed in 
Table|5] Using these four short GRB opening angle measure¬ 


ments alone, taking into account the published range of angles 
for individual bursts, the median is {6j) =6zt I deg (Figure[ 8 ]l. 


Table 5 

Short GRB Opening Angles 


GRB 

Band" 

Oj 

(deg) 

5tL 

(days) 

Reference 

050709 

0 

> 15° 

16.2 

1 

050724A 

X 

>25° 

22.0 

2 

051221A 

X 

6-7° 

26.6 

3 

090426A 

0 

5-7° 

2.7 

4 

101219A 

X 

> 4 ° 

3.9 

5, This work 

111020A 

X 

3-8° 

10.2 

6 

111117A 

X 

> 3-10° 

3.0 

7, 8 

120804A 

X 

> 13° 

45.9 

9, This work 

130603B 

OR 

4-8° 

6.5 

10 

140903A 

X 

^ 6 ° 

3.0 

11, This work 

140930B 

X 

> 9 ° 

23.1 

This work 


Note. — Bursts with opening angle measurements are in bold. 

“ This indicates the band in which the jet break was detected or the 
lower limit was placed. X=X-ray, 0=optical, R=radio. A range of 
angles is due to uncertainty in the redshift, kinetic energy or circum- 
burst density. 

* This is the time elapsed between burst detection and the last ob¬ 
servation. 


References: (1) iFox et all 120051; 12) IGruoe et al.| 

20061: 13) 

Soderbers et al.l 120061: 14) INicuesa Guelbenzu et al.l 

20111: 15i 

Fons et al.ll2013l: 16)IFone et al.ll2012bl: 17)IMareutti et alJ2012l: 18) 

Sakamoto et al.ll2013l: 19)[Berser et al.ll2013l: 110)lFons 

et alJl20l4 


(lli lEvans et ai.ll 20 (iS 


However, the majority of short GRBs do not have detected 
jet breaks and instead exhibit a single power-law decline as 
long as they are detected. In these cases, the time of the 
last observation (5fiast) can be used to place lower limits on 
the opening angles. The inclusion of these bursts is essen¬ 
tial in understanding the true opening angle distribution. In 
most cases, Swift/XRT observes short GRBs until they fade 
below the detection threshhold at < 1 d ay, and enables rela¬ 
tively shall ow lower limits of > 2-5° (iCoward et al.ll20 1^ 
iFong et ^I2012bh . The inclusion of such limits will not have 
a significant effect on the opening angle distribution, as they 
virtually span the entire range of allowable angles. 

Therefore, in order to have a more complete understanding 
of the opening angle distribution, we collect all existing pub¬ 
lished lower limits of 9j >5°, and calculate lower limits for 
GRBs 101219A, 120804A, 140903A, and 140903B using the 
observations and physical parameters presented in this work. 
The inferred lower limits, the band in which the jet break was 
detected, and the value of (jfiast used to compute the limits, are 
listed in Table |5] These seven events demonstrate that multi¬ 
wavelength afterglow observations to <5fiast ~ 3 - 25 days en¬ 
able more meaningful lower limits on the opening angles of 
>5-25° (Table©. 

To calculate the opening angle distribution, we give each 
of the 11 events equal weighting, where measurements are 
given Gaussian probability distributions to represent their al¬ 
lowed range of angles (Table©, while lower limits are given 
probability that is evenly distributed between the lower limit 
and the maximum possible opening angle, max = 90°. The 
resulting cumulative distribution for 11 short GRBs, includ¬ 
ing measurements and lower limits, is shown in Figure © 
Assuming an upper bound of max = 90°, the short GRB 
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Figure 8. Cumulative histograms of 11 short GRBs (dark red dash-dotted) 
and 265 long GRBs (dark blue) with opening angle measurements or limits, 
assuming a maximum opening angle of 0j^tnax = 90°. Also shown is the distri¬ 
bution for 11 short GRBs assuming a more realistic maximum opening angle 
of ^j,max = 30° (red). The distributions for the subset of four short GRBs 
(light red dotted) and 248 long GRBs (light blue) with opening angle mea¬ 
surements are also shown. Medians are denoted by color-coded an'ows, and 
listed with their Icr uncertainties in units of degrees. 


population median is { 0 j) = deg (Icr). Motivated by 
simulations of pos t-merger black hole accretion predict jets 
with 0,^5- 3 ( )° (iRuffert & .Tankalfp^ I Alov et al.ll200l 
iRosswo jl2005t iRezzolla et al.ll201 if we also calculate the 
cumulative distribution employing a more realistic maximum 
value of 0j max = 30°, and find a median of { 6 j) = 16 ± 10 deg 
(la). 

To compare these distributions to those for long GRBs, 
we collect opening angle measuremen ts for 265 long 


Berger et al.ll2003l iBloom et alJl2003l iGhirlanda 

et al.l 

2004; 

Friedman & BloomI 120051: 

Racusin et al.l 120091: 

Cenko et al. 

20101 2011: Filgasetal 

120111 iGoldstein et al.1 120111; 


Rvan et al.ll2015h and calculate the cumulative distributions 


in the same manner (Figure [ 8 ]l. We find a median value for 
the 248 long GRBs with measurements of (0j) = 13^| deg. 
Including the 17 events with limits (0j_niax = 90°), the median 
becomes 14;!;jg deg. 

The opening angle distribution of short GRBs impacts 
the true energy scale, as the true energy is lower than the 
isotropic-equivalent value by the beaming factor, fh, where 
fh = 1 -cos(0j) and therefore Etme = fhEi^o- To calculate the 
cumulative beaming factor distribution, we use the individ¬ 
ual opening angle probability distributions for each burst to 
convert to individual distributions in beaming factor. We then 
sum the individual distributions in a cumulative sense and cal¬ 
culate the median and Icr uncertainties about the median. In¬ 
cluding all short GRBs with opening angle measurements and 
limits and assuming the more realistic scenario of ^j.max = 30°, 
the median beaming factor is fi, = O.OT^ggj. The beam¬ 
ing correction is less substantial if we assume ^j.max = 90°, 
fb = 0 . 1 ?;!;®-jg, and is much more substantial if we only include 
short GRBs with measurements, fh = 0.005 ± 0.002. 

We find median isotropic-equivalent 7 -ray and kinetic en¬ 
ergy scales of ~ 2 x 10^' erg and Tsk^so ~ (1 - 3) x 
10^' erg. Applying the beaming correction for the most re¬ 


alistic scenario gives median beaming-corrected 7 -ray and ki¬ 
netic energy scales of {Ef) = 0.8;!;Qg x 10^° erg and {Ek) = 
0 . 8 ;!;q 7 X 10^° erg, resulting in a total beaming-corrected en¬ 
ergy release of (istot) = l- 6 ^j 3 x 10^° erg. The inferred energy 
scales can be used to constrain the mechanism of energy ex¬ 
traction to power the relativistic jet: the therma l energy re¬ 
lease from vv annihilation in a b aryonic outflow (iJaroszvnskil 
1 19931 iMochkovitch et all 1 19931) and magnetohydrodynamic 
(MHD) processes in the black hole’s accretion remnant (e.g. 
iBlandford & Znajekll977HRosswog et alJ2003h . The general 
consensus is that vv annihilation can only produce beaming- 
corrected total energy releases of 10"^^ - 10"^^ erg, while MHD 
processes can more easily produce energy releases in excess 
of 10^^ erg (IRuf fert & Jankal 1 199 9Blal: iPoph am et al .j 1 19(^ 

[RosswodioOsHBirkl et alJ2007HLee & Ramirez-Ruizl2007h . 

Thus, if the majority of short GRBs have wider opening an¬ 
gles than the four short GRBs with measurements, and thus 
have a smaller overall correction to the isotropic-equivalent 
energy scale, it will be necessary to invoke MHD processes to 
explain the observed energy releases. 

The opening angles also impact the event rate, as the 
true event rate is elevated compared to the observed rate 
by a factor of /^', so Ktme = The current es¬ 

timated obs erved short GRB volumetric rate is JJobs ~ 10 
Gpc“^ yr“^ (iNakar et al.ll2006h . Using = 27^}g*, which 
corresponds to all short GRBs with opening angle measure¬ 
ments and limits ( 6 *j,niax = 30°), we find a true event rate of 
JRtme ~ 270^}gg° Gpc“^ yr“^. The observed all-sky event rate 
of « 0.3 yr“^ within 200 Mpc (iGuetta & PiranI l2005h then 
becomes 8 j^^ yr“'. We note that this rate is conservative 
compared to previously reported rates based on short GRB 
observations dCoward et alJ 120121 : iFong et ^ 1201 2bl 120141) . 
as it properly incorporates opening angle lower limits, with 
the only assumption of an upper bound on the opening an¬ 
gle of 30°. This range is fully consistent with the expected 
detection rates of neutron star mergers within a volume of 
200 Mpc by Advanced LIGOWIRGO of Ri 0.2 - 200 yr'^ 
(!LIG0 Scientific Collaboration et al.ll2013l) . Since there are 
a limited number of short GRBs with meaningful informa¬ 
tion on the opening angles, any additional measurements will 
greatly help to elucidate the true opening angle distribution, 
and therefore true energy scale and event rate. 

5.3. Connection to Galactic Environments 

We connect the afterglow properties of short GRBs to their 
larger-scale, galactic environments. In particular, we inves¬ 
tigate trends between the circumburst densities, which probe 
the explosion environment on sub-parsec-scales, to the pre¬ 
dicted distributions for NS-NS mergers, host galaxy type, and 
locations within the host galaxies. 

We find a wide range of inferred circumburst densities 
for the 38 bursts that we have studied in detail, and can 
compare the properties of the sub-parsec-scale environment 
with the global host galaxy environmen t. Separat i ng the 
bursts by host galaxy type according to IFong et aH (120131) 
and additional data collected since, we find that bursts in 
both star-forming and elliptical host galaxies span a wide 
range of densities, 10“®-5 cm“^ (Figure |9]). We compare 
the inferred densities to predictions for NS-NS mergers from 
population synthesis for varying Galactic potentials, which 
have in put distributions for me rger timescales and kick ve¬ 
locities (iBelczvnski et al.l[2006h . Considering only the four 





































































Short GRB Afterglows 


19 


1- 

* -1 

^ mu —!—1 1 )K 1 

l— ^ 

•s 

-j-H «H 1*1 

r — ' 1 — ^ i-#j 

N 

T"* SF* 

^ Ell ■ 


10 -^ 



Figure 9. Top: Projected physical offset, 5R, versus circumburst density 
for 22 short GRBs with inferred circumburst densities and sub-arcsecond 
localization allowing for precise offset measurements. The sample is sep¬ 
arated by host galaxy type: star-forming hosts (blue asterisks) a nd elliptical 
hosts with no signs of star formation (red squares) according to lFong et alj 
For bursts with no spectroscopic redshift, we assume z = 0.5 to cal¬ 
culate the physical offset, except in two cases, GRBs 060313 and I11020A, 
where we have assumed z = 1 as suggested by the afterglow modeling. 
Middle: Host-normalized offset, SRjre, versus circumburst density for 16 
bursts with measured host gala x y sizes. Physical and host-normalized off ¬ 
sets are fromIFong et alJ i20T(il). iMargutti et alJ 1201211 . iBerger et af] 1201311 . 
IFong & Bergej 1201311 . ISakamotoer^jl20T3l) 7and this work. En'or bars cor¬ 
respond to 1(T confidence. Bottom: Histogram of inferred median circum¬ 
burst densities for 22 short GRBs with physical offsets. 

bursts with elliptical hosts, we find particularly good agree¬ 
ment with the distributions for the large elliptical galaxy 
model (M* = 1 O"M 0 , Mhaio = IO'^Mq) which has probability 
peaks at 10 “^ cm“^ and 1 cm“^, and poor agreement with the 
small elliptical galaxy model (M* = IO^Mq), which is domi¬ 
nated by very low densities of < 10“® cm“^ (iBelczvnski et al.l 
120061) . This conclusion is commensurate with the known stel¬ 
lar masses of short GRB elliptical hosts, which have a median 
of M* R:! 10"" Mq (iLeibler&BergeillMol: iBergeillMl . 
There are no predicted density distributions for star-forming 
galaxies; however, we would expect short GRBs which origi¬ 
nate in elliptical galaxies to have lower inferred densities due 
to the lower average I SM densities in elliptical galaxies (e.g., 
iFukazawa et akllTOObl) . Based on the small number of events, 
we do not find any significant difference between the inferred 
circumburst densities of short GRBs in star-forming versus 
elliptical hosts. 

We next investigate the relationship between the inferred 
circumburst densities and burst offsets from their host galax¬ 
ies. If short GRBs trace the large-scale distribution of the 
ISM, we expect the inferred circumburst densities to decrease 
as a function of offset. In this vein, we gather all available pro¬ 
jected physical offsets, 5R, between the aftergl ow location and 
host galaxy center, derived from ground-based (iMargutti et al.l 


120121: IBerger et al.l 120131: ISakamoto et al.l 120131) and Hub- 
ble Space Telescope (HST) observations (iFong et al.l 120101: 
iFong & Bergerll2013l) . For bursts with no spectroscopic red- 
shift, we calculate the physical offsets at z = 0.5 to be con¬ 
sistent with this work. However, for GRBs 060313 and 
111 020A, afterglow modeling implies that z > 0 .5 (c.f. Sec¬ 
tion [JjJ IFong et ani2012bl: iRoming et al.ll200^ so we keep 
the original fiducial value of z = 1 for these bursts for complete 
uniformity. The distribution of circumburst densities with re¬ 
spect to projected physical offset for 22 bursts is shown in Fig- 
ure|9] We find that four bursts with 6R > 10 kpc have very low 
densities of < 10 “^ cm“^, while bursts with < 10 kpc have 
a wider range of densities, spanning Ri 6 orders of magnitude 
(Figure |9]l. Overall, we find that for 6R Ki 1-15 kpc, there is 
no obvious trend between circumburst density and projected 
physical offset. We also find no obvious trends when con¬ 
sidering only bursts with relatively well-measured densities 
(Vc < Vx)- 

To analyze the relationship with offsets in a more uniform 
manner, we utilize offsets that have been normalized by the 
sizes of their host galaxies (5R /where is the galaxy half- 
light radius). The sample of bursts with host-normalized off¬ 
sets is smaller since precise galaxy size measurements require 
the resolution of HST, thus the sample c omprises 16 events 
(IFong et al.l 120101: IFong & Berged 12013b . The circumburst 
densities as a function of projected host-normalized offset is 
provided in Figure|9] Our analysis suggests that for < 5 r^, the 
inferred densities are largely independent of host-normalized 
offset. We discuss a couple of possible contributing factors. 
First, since we can only measure projected offsets, we are not 
sensitive to the distance component along our line of sight, 
which could contribute a significant amount to the absolute 
distance. This may explain the case of GRB 061006, which 
has a small projected offset Ri QAre but has a very low density 
of R:^ 2 X 10“^ cm“^ (Table [3] and Figure |9]l. Second, the af¬ 
terglow only probes the sub-parsec circumburst environment 
and to a certain extent will be more sensitive to small-scale 
fluctuations in the ISM rather than the average ISM density 
on kiloparsec scales. 

However, bursts that appear to have no coincident host 
galaxy to deep opticaFNIR limits of > 26 mag and are lo¬ 
cated Ri 30-75 kpc fro m the nearest most probable host 
galaxy (“host-less” bup ts ; lBergeri 2 ()Tol: IFong & Berged2013l: 
iTunnicliffe et al.l2013b are expected to have low inferred den¬ 
sities. Indeed, the three bursts located at offsets of > lOr^ 
have low densities of < 10 “^ cm“^, as expected if these bursts 
occur in the IGM or outer halos of their hosts. 

6. CONCLUSIONS 

We present the most comprehensive catalog of short GRB 
afterglows to date, representing a decade of observations since 
the launch of Swift in 2004. This catalog is comprised of 103 
short GRBs with prompt X-ray, optical/NIR and radio follow¬ 
up, enabled by broad-band Target-of-Opportunity programs. 
Applying the synchrotron afterglow model to the observa¬ 
tions, we also place constraints on the isotropic-equivalent 
kinetic energies and circumburst densities for a subset of 38 
events with well-sampled data sets. While a handful of short 
GRB afterglows have been studied in detail on an individual 
basis, our work presents the energy and density scales for a 
large population of events for the first time. We come to the 
following key conclusions: 

• The afterglow observations presented in this work in- 
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elude 71 X-ray detections, 30 optical/NIR detections, 
and 4 radio detections. The detection fractions are 91%, 
40% and 7%, respectively, after accounting for observ¬ 
ing constraints. We present new optical/NIR observa¬ 
tions for 11 events, and new radio observations for 25 
bursts. 

• Applying a synchrotron model to the broad-band after¬ 
glows of 38 bursts, we calculate the inferred circum- 
burst density. Considering a range of scenarios with 
varying values for the microphysical parameters, cool¬ 
ing frequency, and minimum circumburst density, the 
median circumburst density is (3 - 15) x 10“^ cm“^, 
with 90% upper limits of ugo ~ 0.4-3 cm“^. Further¬ 
more, « 80-95% of the probability is below Ri 1 cm“^. 
This indicates that overall short GRBs explode in low- 
density environments. 

• Depending on the set of assumptions in our analy¬ 
sis, we infer a median isotropic-equivalent kinetic en¬ 
ergy of « (1-3) X 10^' erg (considering all scenar¬ 
ios), and an isotropic-equivalent 7 -ray energy scale of 
Ri 2 X 10^' erg. We find a median 7 -ray efficiency of 
Ri 0.40-0.56. 

• We find no obvious trends between circumburst density 
and host galaxy offset for projected offsets of < 10 kpc 
(or < 5 re), and no trend between density and host 
galaxy type, indicating that the circumburst density is 
not strongly dependent on the average ISM density. 
However, three bursts in our sample with offsets of 
>10 kpc have low densities of < 10 ““* cm“^, as ex¬ 
pected if these bursts explode in the IGM. 

• Using 11 short GRBs with opening angle measure¬ 
ments and lower limits, and assuming a maximum value 
on the opening angle of 30°, we calculate a median jet 
opening angle of 16 ± 10 deg and a median beaming 
factor of 0.04 ^q Q 3 . This results in a beaming-corrected 
total energy release of Ri ^ x lO^** erg (Icr range), 
which is broadly consistent with the two primary pro¬ 
posed mechanisms of energy extraction, vv annihila¬ 
tion and MHD processes. The beaming-corrected vol¬ 
umetric rate is Ri 270^}gg° Gpc“^ yr“* with an all-sky 
event rate within a volume of 200 Mpe of 8 ^ 15 ^ yr“'. 
This range is fully consistent with the expected detec¬ 
tion rates of gravitational wave signals from neutron 
star mergers by Advanced LIGO/VIRGO within the 
same volume. 

Our study highlights the importance of broad-band obser¬ 
vations in constraining the basic properties of short GRBs. 
For bursts with detected radio afterglows, we can start to con¬ 
strain the microphysical parameters, which has thus far only 
been possible for long GRBs. While our study provides the 
isotropic-equivalent 7 -ray and kinetic energy scales, the true 
energy release depends on the degree of jet collimation for 
short GRBs. Current knowledge of the collimation of short 
GRBs relies on only a handful of events with measured open¬ 
ing angles from their light curves, primarily due to the faint¬ 
ness of short GRB afterglows which prevent temporal mon¬ 
itoring on timescales longer than 1-2 days. Therefore, it 
is imperative to use the most sensitive ground- and space- 
based resources to uncover additional collimated events or 


place meaningful lower limits on the opening angles. It is 
especially important to undertake these studies while Swift is 
in operation, since this satellite has the unique capability of 
providing multi-wavelength light curves within minutes after 
the bursts. 

The past decade of short GRB observations has enabled 
significant progress in understanding the basic properties 
of short GRBs, namely their energetics, circumburst den¬ 
sities, and opening angles. Furthermore, in addition to 
informing the behavior of on-axis afterglows, the circum¬ 
burst density and energy are key parameters which feed 
in to predictions for electromagnetic counterp arts to com¬ 
pact object mergers, such as off-axis afterglows (iGranot et alJ 
l2002t Ivan Eerten et al.ll2010l) and long-lived rad io flares from 
mildly relativistic ejecta (iNakar & PiranI 120111) . Advanced 
LIGO/VIRGO will dete ct NS-NS mergers within a horizon 
distan ce of rs 200 Mpe (ILIGO Scientific Collaboration et al.l 
120131) . making these alternative electromagnetic signatures 
promising for joint detection with gravitational waves. In 
a subsequent work, we will assess the detectability of such 
counterparts by using the distributions of circumburst densi¬ 
ties and energies of on-axis short GRBs as inputs, which will 
help to inform search strategies in the upcoming revolutionary 
era of gravitational wave discovery. 
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APPENDIX 

BROAD-BAND AFTERGLOW CATALOGS 

Table Al 

Short GRB Late-time X-ray afterglow catalog 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(MJy) 

o-x 

(/iJy) 

References 

050509B 

831 

381 

0.07 

0.02 

Gehrels et al. (2005) 


6.2 X 10^ 

2.0 X 10* 

<0.02 



+2.1 X 10^ 

5.0 X 10'* 

< 8.0 X 10^* 



050709 

1.4 X 10^ 

1.5 X lO'* 

1.9 X 10“* 

1.4 X 10^ 

Fox et al. 120051 


2.1 X 10^ 

5.0 X 10* 

< 1.5 X 10^* 



+2.2 X 10^ 

3.8 X lO'* 

1.1 X 10“* 

2.0 X 10'* 



2.8 X 10^ 

6.4 X 10* 

< 7.3 X 10-+^ 




3.7 X 10^ 

2.1 X 10* 

< 5.9 X 10“* 




+ 1.4 X 10* 

6.1 X 10* 

1.3 X 10-* 

6.0 X 10'* 



+ 1.4 X 10* 

1.2 X lO"* 

1.1 X 10-^ 

2.0 X 10'* 


050724A 

1.0 X 10^ 

514 

0.43 

0.086 

Gmoe et al. (2006) 


1.4 X 10^ 

323 

0.68 

0.14 


1.9 X 10^ 

493 

0.44 

0.090 



6.6 X 10^ 

2.0 X 10* 

0.091 

0.022 



1.3 X 10“* 

3.4 X 10* 

0.050 

0.013 



2.5 X 10“* 

1.9 X 10* 

0.10 

0.023 



3.5 X 10“* 

805 

0.27 

0.056 



3.7 X lO'* 

1.2 X 10* 

0.19 

0.039 



4.1 X lO’* 

830 

0.26 

0.053 



4.2 X 10"* 

652 

0.32 

0.068 



4.2 X lO’* 

772 

0.27 

0.057 



4.7 X 10"* 

1.2 X 10* 

0.18 

0.038 



4.8 X 10^ 

880 

0.24 

0.050 



4.8 X 10^ 

790 

0.27 

0.056 



5.2 X 10"^ 

587 

0.37 

0.075 



5.3 X lO’* 

800 

0.27 

0.055 



5.8 X 10^ 

759 

0.29 

0.058 



5.9 X lO’* 

554 

0.39 

0.080 



6.3 X 10^ 

972 

0.22 

0.046 



6.6 X 10“* 

1.0 X 10* 

0.20 

0.044 



7.7 X lO’* 

1.6 X 10* 

0.11 

0.027 



8.8 X 10^ 

2.0 X 10* 

0.10 

0.023 



9.9 X 10"* 

1.2 X 10* 

0.090 

0.028 



1.4 X 10* 

4.4 X 10* 

0.025 

8.4 X 10'* 



1.5 X 10* 

4.0 X 10* 

0.027 

9.5 X 10'* 



1.8 X 10* 

7.0 X 10* 

0.012 

5.5 X 10'* 



2.2 X 10* 

2.1 X 10* 

0.014 

8.4 X 10'* 



5.6 X 10* 

4.7 X lO'* 

1.4 X 10-* 

7.8 X 10^ 



+2.1 X 10* 

8.3 X 10* 

6.0 X 10'* 

1.1 X 10'* 



+2.1 X 10* 

9.6 X 10* 

5.1 X 10'* 

9.7 X 10^ 



+2.3 X 10* 

1.2 X lO'* 

3.9 X 10'* 

7.5 X 10-+^ 



+2.4 X 10* 

1.3 X lO'* 

3.8 X 10'* 

7.2 X 10-+^ 



+2.5 X 10* 

6.5 X 10* 

5.4 X 10'* 

1.2 X 10'* 



+ 1.9 X 10* 

4.3 X 10-* 

3.3 X 10'^ 

1.2 X 10-+^ 


050813 

4.7 X 10* 

2.7 X 10* 

3.4 xlO'* 

1.0 xlO'* 

Evans et al. (2007a. 20091 










24 


Fong et al. 


Table Al — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(^Jy) 

o-x 

(/^Jy) 

References 


6.5 X 10^ 

3.4 xlO* 

<2.7 X 10^3 



050906 

3.0 xlO^ 

5.2 xlO^ 

<0.018 


Pagamet^ Q0051 

050925 

1.4 xlO^ 

1.4 xlO^ 

<6.6 X 10'3 


Beardmore et al. (20051 

051210 

894 

840 

0.40 

0.15 

This work 

051221A 

1.4 X 10^ 

512 

0.86 

0.10 

Soderbers et al. f2006k Belczvnski et al. (2006) 


1.9 X 10^ 

512 

0.59 

0.085 



2.4 X 10^ 

512 

0.43 

0.082 



5.3 X 10^ 

3.5 X lO^* 

0.23 

0.033 



8.8 X 10^ 

3.5 X 10^ 

0.24 

0.035 



1.2 X 10“* 

3.5 X 103 

0.25 

0.026 



1.6 X 10^ 

3.5 X 103 

0.20 

0.076 



1.9 X 10“* 

3.5 X 103 

0.16 

0.021 



2.3 X 10^ 

3.5 X 103 

0.11 

0.023 



2.6 X 10^ 

3.5 X 103 

0.11 

0.026 



3.0 X 10^ 

3.5 X 103 

0.088 

0.015 



3.7 X 10“* 

3.5 X 103 

0.13 

0.020 



4.5 X 10“* 

1.0 X lO'* 

0.079 

0.011 



5.5 X 10^ 

1.0 X lO”* 

0.066 

0.010 



1.0 X 10^ 

4.1 X lO'* 

0.017 

2.7 X 10-3 



1.4 X 10^ 

4.1 X 10“' 

0.018 

2.9 X 10-3 



2.0 X 10^ 

5.7 X 10“* 

0.013 

1.9 X 10-3 



2.5 X 10^ 

5.7 X 10“' 

0.014 

2.0 X 10-3 



3.1 X 10^ 

5.7 X 10“' 

9.0 X 10^3 

1.7 X 10-3 



3.8 X 10^ 

8.6 X 10“' 

5.2 X 10-3 

1.0 X 10-3 



4.7 X 10^ 

8.6 X 10“' 

5.8 X 10-3 

1.1 X 10-3 



5.5 X 10^ 

8.6 X 10“' 

4.4 X 10-3 

9.4 X 10-+' 



6.4 X 10^ 

8.3 X 10“' 

3.6 X 10-3 

9.4 X 10-+' 



7.3 X 10^ 

8.2 X 10“' 

5.4 X 10-3 

1.0 X 10-3 



9.8 X 10^ 

4.1 X 10* 

2.3 X 10-3 

4.1 X 10-+' 



+ 1.3 X 10^ 

5.4 X 103 

0.022 

2.2 X 10-3 



+ 1.4 X 10^ 

5.5 X 103 

0.021 

2.1 X 10-3 



+ 1.4 X 10^ 

6.4 X 103 

0.018 

1.9 X 10-3 



+ 1.5 X 10^ 

6.6 X 103 

0.018 

1.8 X 10-3 



+ 1.6 X 10^ 

6.3 X 10' 

0.017 

1.8 X 10-3 



+3.9 X 10^ 

9.5 X 103 

6.1 X 10-3 

8.8 X 10-+' 



+4.0 X 10^ 

1.0 X 10“' 

5.5 X 10-3 

8.0 X 10-+' 



+4.1 X 10^ 

1.0 X 10“' 

5.3 X 10-3 

7.9 X 10-+' 



+ 1.3 X 10* 

1.8 X 10“' 

5.5 X 10-“' 

1.9 X 10-+' 



+ 1.7 X 10* 

2.5 X 10“' 

3.8 X 10-“' 

1.4 X 10-+' 



+2.3 X 10* 

4.8 X 10“' 

1.9 X 10-“' 

7.2 X 10-3 



060121 1.1 x 10“' 

153 

0.61 

0.14 

Evans et al. (2007a. 20091 

1.1 X 10“' 

221 

0.42 

0.095 


1.1 X 10“' 

206 

0.45 

0.10 


1.1 X 10“' 

198 

0.47 

0.11 


1.1 X 10“* 

208 

0.34 

0.087 


1.2 X 10“' 

178 

0.52 

0.12 


1.2 X 10“* 

183 

0.51 

0.12 


1.2 X 10“' 

155 

0.60 

0.14 


1.2 X 10“' 

193 

0.36 

0.094 


1.2 X 10“' 

338 

0.37 

0.072 


1.6 X 10“' 

281 

0.27 

0.069 


1.6 X 10“' 

308 

0.24 

0.063 


1.6 X 10“' 

196 

0.39 

0.10 


1.7 X 10“* 

258 

0.29 

0.076 


1.7 X 10“* 

246 

0.31 

0.080 


1.7 X 10“' 

236 

0.32 

0.083 


1.8 X 10“' 

221 

0.34 

0.089 


1.8 X 10“' 

201 

0.38 

0.10 


1.8 X 10“' 

617 

0.19 

0.039 


2.2 X 10“' 

361 

0.23 

0.056 


7.0 X 10“' 

1.3 X 103 

0.060 

0.015 


7.5 X 10“' 

2.5 X 103 

0.030 

7.9 X 10-3 


8.1 X 10“* 

2.3 X 103 

0.053 

0.011 


8.9 X 10“' 

7.6 X 103 

0.043 

8.0 X 10-3 


1.6 X 10* 

1.2 X 10“' 

0.021 

4.9 X 10-3 
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Table A1 — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(^Jy) 

o-x 

(/^Jy) 

References 


1.7 X 10^ 

1.9 X 10® 

0.038 

0.010 



1.8 X 10^ 

7.3 X 10® 

0.033 

7.2 X 10'® 



2.0 X 10^ 

3.6 X 10® 

0.012 

3.2 X 10'® 



2.5 X 10^ 

4.9 X 10® 

0.013 

3.4 X 10'® 



3.1 X 10^ 

6.5 X 10® 

7.7 X 10^® 

2.1 X 10'® 



4.2 X 10^ 

1.5 X 10® 

4.1 X 10“® 

1.2 X 10'® 



6.3 X 10^ 

2.8 X 10® 

3.5 X 10^® 

9.7 X 10-^ 



1.0 X 10® 

8.3 X 10® 

1.8 X 10“® 

8.5 X 10-^ 


060313 

4.2 X 10^ 

161 

1.3 

0.52 

This work 


4.3 X 10^ 

92 

1.4 

0.52 



4.4 X 10^ 

97 

1.7 

0.62 



4.5 X 10^ 

128 

1.2 

0.43 



4.7 X 10^ 

113 

1.2 

0.45 



4.8 X 10^ 

113 

1.1 

0.39 



4.9 X 10^ 

128 

1.1 

0.43 



5.0 X 10^ 

65 

2.4 

0.91 



5.1 X 10^ 

133 

0.91 

0.33 



5.2 X 10^ 

146 

0.83 

0.30 



5.4 X 10^ 

176 

0.64 

0.23 



5.5 X 10^ 

163 

0.75 

0.27 



5.7 X 10^ 

95 

1.4 

0.54 



5.8 X 10^ 

211 

0.61 

0.23 



6.0 X 10^ 

125 

1.2 

0.46 



6.1 X 10^ 

153 

0.85 

0.31 



6.3 X 10^ 

125 

1.0 

0.38 



6.5 X 10^ 

286 

0.96 

0.30 



1.0 X lO'* 

301 

0.38 

0.14 



1.0 X 10“* 

309 

0.43 

0.16 



1.1 X 10"* 

228 

0.66 

0.25 



1.1 X 10“* 

243 

0.47 

0.17 



1.1 X 10“* 

266 

0.44 

0.16 



1.1 X lO'* 

384 

0.346 

0.13 



1.2 X 10“* 

344 

0.36 

0.13 



1.2 X lO'^ 

459 

0.37 

0.12 



1.6 X 10“* 

602 

0.19 

0.069 



1.6 X 10^ 

412 

0.29 

0.11 



1.7 X 10“* 

414 

0.27 

0.10 



1.8 X 10^ 

1.1 X 10® 

0.25 

0.078 



2.2 X 10"* 

934 

0.11 

0.043 



2.3 X 10“* 

545 

0.20 

0.075 



2.3 X 10^ 

1.1 X 10® 

0.13 

0.047 



2.8 X 10“* 

781 

0.15 

0.057 



2.9 X lO'* 

1.8 X 10® 

0.089 

0.031 



3.4 X 10“* 

2.5 X 10® 

0.061 

0.022 



4.2 X 10"* 

6.3 X 10® 

0.058 

0.024 



5.2 X 10“* 

1.2 X 10® 

0.036 

0.014 



8.7 X 10^ 

5.9 X 10® 

0.013 

5.4 X 10'® 



2.7 X 10® 

3.1 X 10® 

3.6 X 10-® 

1.2 X 10'® 


060502B 

3.8 xlO'* 

6.6 X10® 

<1.0 xlO-® 


Evans et al. 120073^0091 

060801 

918 

954 

0.20 

0.077 

This work 

061006 

3.6 X 10® 

670 

0.12 

0.067 

This work 


4.4 X 10® 

997 

0.11 

0.060 



9.4 X 10® 

4.9 X 10® 

0.046 

0.026 



1.4 X 10"^ 

4.7 X 10® 

0.054 

0.030 



6.6 X 10“® 

2.4 X 10® 

6.0 X 10^® 

3.3 X 10'® 


061201 

5.7 X 10® 

238 

0.55 

0.31 

This work 


5.9 X 10® 

188 

0.56 

0.31 



6.1 X 10® 

238 

0.45 

0.26 



6.4 X 10® 

264 

0.66 

0.34 



1.2 X 10® 

5.9 X 10® 

0.11 

0.072 



2.1 X 10® 

1.2 X 10® 

0.055 

0.036 



6.7 X 10® 

1.3 X 10® 

4.8 X 10'® 

2.9 X 10'® 


061210 

2.3 X 10® 

5.3 X 10® 

0.048 

0.013 

Evans et al. t2007a. 20091 


3.1 X 10® 

9.4 X 10® 

0.026 

8.0 X 10'® 



4.2 X 10® 

1.3 X 10® 

0.018 

6.0 X 10'® 



7.0 X 10® 

4.9 X 10® 

< 0.014 



061217 

6.0 X10® 

3.2 xlO® 

<2.5 XlO'® 


Evans et al. (2007a. 20091 
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Table Al — Continued 


GRB 

St 

Exposure Time 

Fx 

o-x 

References 


(s) 

(s) 

(^Jy) 

(/iJy) 



070209* 

4.4 X 10^ 

1.9 X 10"* 

< 8.3 X 10^ 


Stratta et al. (2007bl 

070406* 

9.7 X 10“* 

1.3 X 10'* 

< 1.8 X 10“* 


Troja et al. 120071 

070429B 

1.6 X 10^ 

2.0 X10* 

0.32 

0.080 

Evans et al. 12007a. 20091 


1.6 X 10"^ 

3.6 X lO'* 

< 0.024 



OlOlOT 

4.1 xlO"* 

2.2 X lO'* 

0.15 

0.032 

Evans et al. 12007a. 20091 


4.1 xlO^ 

1.9 xlO* 

< 0.032 



070714B 

8.5 X 10^ 

116 

2.0 

0.64 

This work 


2.8 X 10^ 

183 

0.40 

0.18 



3.1 X 10^ 

311 

0.24 

0.11 



3.6 X 10^ 

675 

0.19 

0.09 



5.9 X 10^ 

417 

0.15 

0.075 



6.5 X 10^ 

795 

0.12 

0.058 



9.2 X 10^ 

1.4 X 10* 

0.054 

0.027 



2.9 X 10“* 

5.8 X 10“ 

3.0 X 10“* 

1.7 X 10'* 


070724A 

4.6 X 10^ 

4.8 X 10* 

0.055 

0.032 

This work 


1.1 X 10“* 

7.6 X 10* 

0.026 

0.014 



1.8 X 10“* 

5.2 X 10* 

0.024 

0.013 



2.6 X 10** 

1.2 X I0“ 

0.014 

7.5 X 10'* 



5.3 X lO'^ 

5.5 X 10“ 

0.010 

5.4 X 10'* 



1.8 X 10* 

3.0 X 10* 

1.8 X 10-* 

5.9 X 10^ 


070729 

2.0 X10“* 

6.6 X10“ 

<3.8 X 10'* 


Evans et al. 12007a. 20091 

070809 

1.1 X 10* 

271 

0.27 

0.13 

This work 


1.6 X 10* 

660 

0.20 

0.080 



5.5 X 10* 

321 

0.23 

0.11 



5.9 X 10* 

507 

0.14 

0.068 



6.3 X 10* 

366 

0.23 

0.11 



6.8 X 10* 

555 

0.13 

0.062 



7.4 X 10* 

673 

0.12 

0.058 



1.1 X 10"* 

673 

0.10 

0.047 



1.2 X 10“* 

803 

0.093 

0.044 



1.3 X 10^ 

989 

0.085 

0.039 



1.8 X 10"* 

2.6 X 10* 

0.038 

0.017 



2.4 X 10^ 

2.0 X 10* 

0.050 

0.021 


070810B* 

983 

7.2 X 10* 

< 6.6 X 10-^ 


Starling et al. (20071 

071017* 

2.0 X10* 

2.1 XlO* 

<5.9 XlO'* 


Evans et al. (2007b) 

071IHB" 

3.7 xlO* 

2.2 XlO* 

< 0.023 


Pern et al. (20071 

071227 

6.2 X 10* 

9.5 X 10* 

0.022 

9.8 X 10'* 

This work 


5.7 X 10“* 

2.5 X 10* 

2.6 X 10“* 

8.6 X 10-^ 



080121“ 

2.0 X 10* 

2.2 X 10“ 

<4.0 X 10'* 


Troja et al. (20081 

081023 

1.4 X 10* 

1.7 X 10* 

0.072 

0.024 

This work 


1.8 X 10“ 

5.3 X 10“ 

9.5 X 10'* 

3.2 X 10'* 


080426 

1.1 X 10* 

100 

2.13 

0.48 

Evans et al. (2007a. 20091 


1.2 X 10* 

148 

1.4 

0.32 



1.3 X 10* 

95 

2.2 

0.50 



1.4 X 10* 

108 

2.1 

0.45 



1.5 X 10* 

125 

1.7 

0.38 



1.7 X 10* 

211 

1.01 

0.23 



1.9 X 10* 

191 

1.44 

0.29 



6.2 X 10* 

258 

0.65 

0.17 



6.6 X 10* 

687 

0.24 

0.063 



7.5 X 10* 

978 

0.19 

0.047 



1.3 X 10“ 

2.1 X 10* 

0.11 

0.025 



2.1 X 10“ 

8.1 X 10* 

0.056 

0.012 



3.8 X 10“ 

2.0 X 10“ 

0.021 

5.7 X 10'* 



7.0 X 10“ 

4.2 X 10“ 

4.9 X 10'* 

2.8 X 10'* 


080503 

920 

1.1 X 10* 

0.45 

0.18 

This work 


+''3.9 X 10* 

3.2 X 10“ 

2.1 X 10'* 

9.6 X 10-+^ 

Perley et al. (2009b) 


+ 1.9 X 10* 

3.3 X 10“ 

< 8.0 X 10'* 


Perley et al. (2009bl 

080702A 

3.9 X 10* 

2.2 X 10“ 

0.041 

8.9 X 10'* 

Evans et al. (2007a. 20091 


7.5 X 10“ 

7.6 X 10“ 

<0.013 


080905A 

758 

593 

0.74 

0.31 

This work 

080919 

1.0 XlO* 

1.5 XlO* 

0.39 

0.20 

This work 


3.0 X 10“ 

1.7 XlO* 

5.9 XlO'* 

5.7 X 10'* 


081024A 

597 

1.6 XlO* 

0.39 

0.070 

Evans et al. (2007a. 20091 
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Table A1 — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(MJy) 

o-x 

(/^Jy) 

References 


1.5 xlO'* 

3.7 xlO"* 

<3.2 X 10^3 



081226A 

1.6 xlO”* 

4.3 X 10“* 

0.026 

0.014 

Evans et al. (2007a. 20091 

081226B® 

1.3 xlO"* 

8.0 X103 

<6.6 X 10^3 


Evans & Hoversten (20081 

090305A'’ 

4.8 X 10^ 

1.3 xlO'* 

<2.0 X10-3 


Beardmore et ^^2009) 

090426 

1.0 x 10^ 

95 

1.2 

0.26 

Evans et al. (2007a. 20091 


1.1 X 10^ 

108 

1.0 

0.23 



1.2 X 10^ 

95 

1.2 

0.26 



1.4 X 10^ 

145 

0.77 

0.17 



1.5 X 10^ 

92 

1.2 

0.27 



1.6 X 10^ 

100 

1.1 

0.25 



1.7 X 10^ 

178 

0.62 

0.14 



1.9 X 10^ 

243 

0.45 

0.10 



2.2 X 10^ 

221 

0.68 

0.13 



5.6 X 10^ 

346 

0.27 

0.072 



6.0 X 10^ 

351 

0.27 

0.071 



6.4 X 10^ 

481 

0.20 

0.051 



6.9 X 10^ 

479 

0.20 

0.051 



7.6 X 10^ 

938 

0.18 

0.035 



1.8 X 10“* 

1.2 X lO'* 

0.060 

0.013 



2.3 X 10^ 

7.5 X 10^ 

1.5 X 10-3 

6.5 X 10^ 


090510 

1.0 X 10^ 

31 

2.7 

1.0 

This work 


2.8 X 10^ 

89 

0.61 

0.26 



2.9 X 10^ 

141 

0.47 

0.21 



3.1 X 10^ 

174 

0.37 

0.17 



3.3 X 10^ 

338 

0.18 

0.081 



3.6 X 103 

251 

0.22 

0.10 



3.9 X 103 

321 

0.18 

0.084 



6.2 X 103 

770 

0.051 

0.027 



6.8 X 103 

547 

0.073 

0.039 



1.3 X 10^ 

9.2 X 103 

0.020 

0.010 



4.1 X 10“* 

7.6 X lO'* 

4.0 X 10-3 

1.3 X 10-3 


090515 

8.1 X 10^ 

2.2 X 103 

0.24 

0.093 

This work 

090607 

9.3 X 10^ 

2.0 X I03 

0.060 

0.017 

Evans et al. (2007a. 20091 


2.0 X 10“* 

2.4 X lO'* 

< 6.6 X 10-3 



09062IB 

7.2 X 103 

8.3 X I03 

0.04 

0.01 

Evans et al. (2007a. 20091 


2.5 X 10^ 

2.5 X lO'* 

6.0 xlO-3 

3.4 X 10-3 


090916* 

6.5 X 10'* 

1.2 xlO'* 

< 9.4 X 10^ 


^ojae^^ (20091 

091109B 

5.7 X 103 

2.4 X 103 

0.097 

0.019 

Evans et al. (2007a. 20091 


1.6 X 10“* 

1.3 X lO'* 

0.030 

0.023 



5.8 X 10“* 

5.8 X lO'* 

0.011 

7.9 X 10-3 



2.0 X 10^ 

1.2 X 10^ 

0.011 

9.6 X 10-3 



100117A 

3.1 xlO'* 

3.6 xl03 

2.2 xlO-3 

7.2 xlO^ 

This work 

100206A 

2.4 X 10“* 

9.5 X lO'* 

<3.3 xlO-3 


Evans et al. (2007a. 20091 

100213 

2.0 X lO'* 

3.1 xlO"* 

< 0.033 


Evans et al. (2007a. 20091 

100625A 

1.6 X lO"* 

1.3 X 10“* 

4.4 X 10-3 

1.8 X 10-3 

Fong_et^ (20131 


9.0 X 10^ 

4.6 X lO'* 

< 8.2 X 10-3 


Eong_et^ Q0131 

100628A“ 

+3.9 xl03 

2.0 xlO'* 

<5.9 X 10-* 


This work 

100702A 

1.4 xlO”* 

4.3 X 10“* 

<0.011 


Evans et al. (2007a. 20091 

101219A 

7.7 x 103 

317 

0.58 

0.20 

Fong_et^ QOl^l 


8.6 X 103 

1.4 X KF* 

0.038 

0.014 



+3.4 X 103 

2.0 X I03 

< 1.9 X 10^ 



101224A 

1.9 xl03 

4.0 xlO'* 

0.015 

5.0 xlO-3 

Evans et al. (2007a. 20091 

110112A 

4.6 X 103 

1.0 xl03 

0.083 

0.028 

Fong_e^^ (20^1 


5.4 X 103 

605 

0.14 

0.047 



6.2 X 103 

941 

0.14 

0.043 



1.1 X lO'^ 

2.5 X 103 

0.064 

0.016 



1.9 X 10“* 

6.6 X 103 

0.025 

8.8 X 10-3 



2.8 X 10“* 

1.3 X lO'* 

0.012 

4.3 X 10-3 



4.5 X lO"* 

2.3 X 10“* 

7.3 X 10-3 

3.1 X 10-3 



1.4 X 103 

2.7 X 103 

2.0 X 10-3 

1.0 X 10-3 


110112B* 

6.1 xlO'* 

1.2 xlO'* 

<2.4 X 10-3 


Littlejohns et ^^20n) 

111020A 

1.1 x 103 

520 

1.0 

0.20 

Eongetal. (2012b) 
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Table Al — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(^Jy) 

o-x 

(/^Jy) 

References 


5.8 X 10^ 

2.5 X 10’ 

0.19 

0.039 



1.2 X 10^ 

2.5 X 10’ 

0.13 

0.032 



1.9 X 10“* 

6.4 X 10’ 

0.11 

0.027 



2.6 X lO'^ 

6.3 X 10’ 

0.14 

0.037 



3.2 X 10^ 

5.9 X 10’ 

0.12 

0.031 



4.2 X 10“* 

1.6 X 10“ 

0.10 

0.028 



1.0 X 10^ 

1.5 X 10’ 

0.019 

6.4 X 10-’ 



2.9 X 10^ 

2.1 X 10’ 

0.013 

4.9 X 10-’ 



++6.9 X lO'^ 

1.4 X 10“ 

0.031 

2.4 X 10-’ 



+2.6 X 10^ 

2.0 X 10“ 

7.0 X 10^’ 

1.0 X 10-’ 



5.1 X 10^ 

2.1 X 10’ 

<0.019 




+8.8 X 10^ 

2.0 X 10“ 

< 1.1 X 10^’ 



111117A 

1.0 X 10^ 

178 

0.55 

0.20 

Mar^utti et al. (20121 


1.2 X 10^ 

211 

0.45 

0.15 



5.0 X 10^ 

881 

0.11 

0.035 



5.5 X 10^ 

1.9 X 10’ 

4.0 X 10-’ 

1.4 X 10-’ 



+2.6 X 10^ 

4.0 X 10“ 

4.2 X 10-“ 

1.6 X 10-+^ 


111121A 

4.2 X 10’ 

85 

2.17 

0.47 

Evans et al. (2007a. 20091 


4.3 X 10’ 

133 

1.3 

0.30 



4.5 X 10’ 

143 

1.4 

0.31 



4.7 X 10’ 

228 

1.4 

0.26 



4.9 X 10’ 

113 

1.6 

0.37 



5.1 X 10’ 

196 

1.1 

0.22 



5.3 X 10’ 

246 

0.93 

0.21 



5.6 X 10’ 

326 

1.0 

0.20 



6.0 X 10’ 

133 

1.3 

0.30 



6.2 X 10’ 

143 

1.2 

0.28 



6.3 X 10’ 

160 

1.1 

0.25 



6.5 X 10’ 

158 

1.1 

0.25 



6.7 X 10’ 

206 

1.1 

0.21 



1.0 X 10^* 

155 

1.2 

0.26 



1.0 X 10^ 

183 

1.0 

0.23 



1.1 X 10“* 

384 

0.34 

0.090 



1.1 X 10^ 

188 

0.94 

0.21 



1.1 X 10“* 

181 

0.74 

0.19 



1.1 X 10^ 

243 

0.54 

0.14 



1.1 X 10“* 

223 

0.60 

0.16 



1.2 X 10^ 

268 

0.50 

0.13 



1.2 X 10"* 

657 

0.37 

0.072 



1.6 X 10^* 

985 

0.24 

0.063 



1.8 X 10^ 

1.6 X 10’ 

0.17 

0.043 



2.2 X 10“* 

667 

0.22 

0.057 



2.3 X 10^ 

1.8 X 10’ 

0.12 

0.026 



4.2 X 10“* 

3.7 X 10“ 

0.045 

8.3 X 10-’ 



1.4 X 10’ 

6.5 X 10“ 

< 0.020 



111222A 

9.7 X 10“* 

7.3 X 10’ 

0.047 

9.5 X 10“’ 

Evans et al. (2007a. 20091 


7.8 X 10’ 

5.2 X 10“ 

< 0.027 




4.0 X 10* 

1.7 X 10’ 

< 0.014 




120305A 

2.3 X 10“ 

9.4 xlO“ 

0.016 

8.0 X10-’ 

Evans et al. (2007a. 20091 

120521A 

6.8 X 10’ 

7.3 X 10’ 

< 0.020 

Evans et al. (2007a, 20091 


120630A 

6.6 X 10’ 

1.2 X 10’ 

0.11 

0.024 

Evans et al. (2007a. 20091 

120804A 

4.3 X 10’ 

203 

1.1 

0.23 

Berser et al. (20131 


4.5 X 10’ 

211 

0.82 

0.20 


4.7 X 10’ 

178 

0.85 

0.22 



4.9 X lO’ 

181 

1.1 

0.26 



5.2 X 10’ 

281 

0.54 

0.14 



5.4 X 10’ 

263 

0.58 

0.15 



5.7 X 10’ 

301 

0.50 

0.13 



6.0 X 10’ 

349 

0.43 

0.12 



6.4 X 10’ 

494 

0.34 

0.090 



1.1 X 10“ 

2.6 X 10’ 

0.29 

0.056 



1.7 X 10“ 

2.5 X 10’ 

0.20 

0.051 



2.3 X 10“ 

2.6 X 10’ 

0.11 

0.021 



3.0 X 10“ 

5.7 X 10’ 

0.096 

0.026 



3.4 X 10“ 

2.6 X 10’ 

0.10 

0.021 
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Table A1 — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(^Jy) 

o-x 

(/^Jy) 

References 


4.3 X 10“* 

8.2 X 10* 

0.049 

0.013 



6.2 X 10“* 

2.3 X lO'* 

0.076 

0.016 



1.4 X 10^ 

9.3 X lO'* 

0.048 

0.014 



3.1 X 10^ 

2.8 X 10* 

6.8 X 10“* 

2.2 x 10^* 



5.0 X 10^ 

2.3 X 10^ 

< 0.0456 




+8.1 X 10^ 

2.0 X 

3.5 X 10“* 

6.2 X 10-^ 



++1.6X 10*^ 

2.5 X lO'* 

3.0 X 10-* 

6.2 X 10^ 



+4.0 X 10* 

5.9 xlO'* 

6.6 xlO-'* 

1.3 xl0‘* 

This work 

1208178* 

6.8 xlO”* 

2.8 xlO* 

<5.9 xlO-* 



121226A 

4.3 X 10^ 

589 

0.58 

0.15 

Evans et al. (2007a. 20091 


4.9 X 10^ 

542 

0.62 

0.16 



5.5 X 10^ 

679 

0.49 

0.13 



7.9 X 10^ 

4.7 X 10* 

0.36 

0.079 



6.2 X 10^ 

8.3 X 10“* 

0.030 

8.6 x 10-* 



1.2 X 10* 

6.2 X 10* 

0.034 

0.027 



130313A“ 

2.1 xlO* 

5.3 xlO'* 

< 0.013 


Evans et al. (2007a. 20091 

130515A 

1.2 X 10“* 

3.0 X 10“* 

<0.014 


Evans et al. (2007a. 20091 

130603B 

4.0 X 10* 

33 

3.8 

1.2 

Fon^ et al. (20141 


4.0 X 10* 

41 

3.3 

1.0 



4.1 X 10* 

75 

1.7 

0.54 



4.1 X 10* 

44 

2.9 

0.95 



4.2 X 10* 

44 

3.0 

0.91 



4.2 X 10* 

36 

3.5 

1.2 



4.3 X 10* 

75 

1.8 

0.53 



4.3 X 10* 

44 

3.0 

1.0 



4.4 X 10* 

54 

2.5 

0.71 



4.4 X 10* 

41 

3.2 

1.1 



4.5 X 10* 

65 

2.0 

0.63 



4.5 X 10* 

31 

4.1 

1.3 



4.6 X 10* 

44 

3.0 

0.99 



4.6 X 10* 

52 

2.4 

0.80 



4.7 X 10* 

65 

2.0 

0.68 



4.7 X 10* 

52 

2.6 

0.74 



4.8 X 10* 

98 

1.6 

0.46 



4.9 X 10* 

57 

2.2 

0.72 



4.9 X 10* 

44 

3.0 

0.98 



5.0 X 10* 

59 

2.3 

0.77 



5.1 X 10* 

83 

1.5 

0.43 



5.1 X 10* 

59 

2.5 

0.71 



5.2 X 10* 

93 

1.4 

0.43 



5.3 X 10* 

96 

1.4 

0.45 



5.4 X 10* 

88 

1.6 

0.49 



5.5 X 10* 

72 

1.8 

0.59 



5.5 X 10* 

41 

3.0 

0.98 



5.6 X 10* 

67 

1.9 

0.59 



5.7 X 10* 

85 

1.5 

0.44 



5.8 X 10* 

78 

1.7 

0.50 



5.8 X 10* 

65 

2.0 

0.59 



5.9 X 10* 

75 

1.6 

0.54 



6.0 X 10* 

67 

1.9 

0.56 



6.0 X 10* 

91 

1.4 

0.42 



6.2 X 10* 

114 

1.1 

0.34 



6.2 X 10* 

65 

2.0 

0.59 



6.3 X 10* 

59 

2.1 

0.59 



6.4 X 10* 

91 

1.4 

0.42 



6.5 X 10* 

122 

1.5 

0.40 



9.8 X 10* 

91 

1.3 

0.36 



9.9 X 10* 

143 

0.83 

0.23 



1.0 X 10“* 

166 

0.71 

0.19 



1.1 X 10“* 

1.5 X 10* 

0.57 

0.16 



1.2 X 10“* 

166 

0.65 

0.19 



1.2 X 10“^ 

226 

0.49 

0.13 



1.2 X 10“^ 

205 

0.72 

0.17 



1.6 X 10“^ 

338 

0.29 

0.076 



1.6 X 10“* 

252 

0.38 

0.11 
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Table Al — Continued 


GRB 

St 

(s) 

Exposure Time 

(s) 

Fx 

(^Jy) 

o-x 

(/^Jy) 

References 


1.7 X lO'^ 

2.0 X 102 

0.21 

0.056 



2.1 X 10^ 

450 

0.19 

0.048 



2.2 X 10“* 

463 

0.19 

0.047 



2.8 X 10^ 

1.1 X 102 

0.077 

0.029 



3.6 X 10^ 

5.9 X 10* 

0.038 

0.011 



4.2 X 10“* 

5.5 X 102 

0.054 

0.014 



4.8 X 10^ 

6.7 X 10* 

0.021 

5.6 X 10-2 



1.5 X 10^ 

2.1 X 102 

4.0 X 10“2 

8.8 X 10-+* 



++2.3 X 10’ 

1.9 X 10'* 

2.4 X 10“2 

3.2 X 10-+* 



++5.6 X 10^ 

3.0 X 10'* 

8.2 X 10-'* 

5.0 X 10^ 


130626A'’ 

8.7 xlO^ 

7.2 xl02 

<1.8 xlO-2 


Pa^e & de Pasauale QO 13") 

130716A 

1.2 X 10"* 

2.5 X 10'* 

0.014 

4.0 X 10-2 

Evans et al. (2007a. 20091 


4.0 X 10“* 

1.3 X 10-* 

< 7.2 X 10^2 



130822A 

3.0 xlO'* 

1.9 xlO"* 

<4.6 X10-2 


Evans et al. (2007a. 20091 

130912A 

9.7 X 10^ 

95 

5.6 

1.3 

Evans et al. (2007a. 20091 


1.1 X 10^ 

103 

2.7 

0.60 



1.2 X 10^ 

100 

1.3 

0.30 



1.3 X 10^ 

105 

1.3 

0.29 



1.4 X 10^ 

95 

1.4 

0.31 



1.5 X 10^ 

180 

1.6 

0.26 



5.2 X 10^ 

780 

0.26 

0.068 



6.4 X 10^ 

1.8 X 102 

0.21 

0.04 



1.2 X 10'* 

2.5 X 10* 

0.066 

0.016 



1.8 X 10'* 

2.5 X 102 

0.057 

0.015 



2.8 X 10'* 

1.4 X 10'* 

0.025 

4.9 X 10-2 



1.1 X 10^ 

5.8 X 10'* 

< 3.9 X 10-2 



131004A 

1.8 X 10^ 

306 

13.2 

3.1 

Evans et al. (2007a. 20091 


4.4 X 10^ 

115 

34.8 

8.1 



5.5 X 102 

115 

33.9 

8.1 



6.9 X 102 

158 

24.8 

5.9 



8.8 X 102 

228 

31.9 

5.5 



1.1 X 10'* 

1.5 X 10* 

0.23 

0.049 



2.4 X 10'* 

1.8 X 10'* 

0.055 

0.015 



5.9 X 10'* 

4.7 X 10'* 

0.050 

0.014 



1.1 X 10^ 

4.1 X 10'* 

< 0.017 



131224A'’ 

1.3 xlO'* 

4.0 X102 

<4.1 xlO-2 


Gomoertz et al. ('2013) 

140129B 

9.9 X 102 

65 

3.01 

0.68 

Evans et al. (2007a. 20091 


1.1 X 102 

45 

4.62 

1.01 



1.1 X 102 

60 

3.24 

0.73 



4.7 X 102 

624 

0.27 

0.052 



6.7 X 102 

366 

0.20 

0.064 


140320A 

1.1 X 102 

1.2 X 10'* 

0.13 

0.032 

Evans et al. (2007a. 20091 


2.7 X lO"* 

1.2 X lO"* 

< 0.062 



140402A*’ 

4.2 X lO-* 

4.5 X 102 

<4.2 X 10-2 


Pagan^^O^) 

140414A'’ 

4.1 X 10'* 

1.7 X 10* 

< 0.027 


D’Avanzo et al. (20141 

140516A 

1.5 X 102 

2.5 X 10'* 

0.013 

2.9 X 10-2 

Evans et al. (2007a. 20(191 


1.9 X 10^ 

3.1 X 102 

1.7 X 10-2 

9.7 X 10^ 


140606A'’ 

4.6 X 102 

2.4 X 102 

< 8.7 X 10-2 


Stroh et al. (20141 

MOeWB* 

5.9 X 10'* 

3.9 X 102 

< 3.4 X 10-2 


Maselli & D’Avanzo (20141 

140622A 

1.3 X 102 

50 

1.86 

0.40 

Evans et al. (2007a. 20091 


2.3 X 102 

3.5 X 10'* 

0.015 

2.9 X 10-2 



8.0 X 10'* 

2.3 X 10'* 

< 8.4 X 10-2 



140903A 

1.0 X 102 

221 

1.02 

0.23 

Evans et al. (2007a. 20091 


1.2 X 102 

176 

1.3 

0.30 



1.4 X 102 

303 

0.74 

0.17 



5.0 X 102 

223 

0.93 

0.16 



5.4 X 102 

158 

0.78 

0.18 



5.6 X 102 

243 

0.95 

0.16 



6.2 X 102 

130 

0.95 

0.21 



6.3 X 102 

123 

0.86 

0.19 



6.5 X 102 

183 

1.05 

0.21 



6.8 X 102 

288 

0.93 

0.16 



7.2 X 102 

338 

0.86 

0.14 



1.1 X 10'* 

160 

0.60 

0.14 



1.1 X 10'* 

216 

0.42 

0.10 
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Table A1 — Continued 


GRB 

St 

Exposure Time 

Fx 

o-x 

References 


(s) 

(s) 

(^Jy) 

(/iJy) 



1.1 X lO'^ 

178 

0.54 

0.12 



1.1 X 10“* 

173 

0.59 

0.13 



1.1 X lO'^ 

216 

0.55 

0.11 



1.3 X lO'* 

333 

0.49 

0.085 



1.6 X lO'* 

223 

0.42 

0.098 



1.7 X 10^ 

228 

0.32 

0.085 



1.7 X 10“* 

203 

0.370 

0.096 



1.7 X 10“* 

133 

0.75 

0.17 



1.7 X 10“* 

333 

0.38 

0.076 



2.2 X 10“* 

517 

0.25 

0.065 



2.3 X 10^ 

354 

0.36 

0.093 



2.3 X lO'^ 

461 

0.27 

0.070 



2.8 X 10“* 

354 

0.23 

0.057 



2.8 X 10“* 

253 

0.30 

0.077 



2.9 X 10“* 

366 

0.20 

0.054 



3.1 X 10“* 

5.1 X 10® 

0.18 

0.033 



4.2 X 10“* 

6.0 X 10® 

0.13 

0.034 



4.6 X lO'* 

479 

0.21 

0.052 



5.4 X 10"* 

6.5 X 10® 

0.080 

0.018 



6.8 X 10^ 

1.2 X lO'* 

0.076 

0.0142 



1.2 X 10^ 

8.6 X lO"* 

0.029 

6.5 X 10'® 



2.6 X 10^ 

5.8 X KF* 

<0.011 



140930B 

1.1 X 10^ 

158 

0.56 

0.13 

Evans et al. (2007a. 20091 


1.2 X 10^ 

115 

0.77 

0.17 



1.3 X 10^ 

133 

0.66 

0.15 



1.5 X 10^ 

196 

0.45 

0.10 



1.7 X 10^ 

191 

0.35 

0.090 



1.9 X 10^ 

221 

0.40 

0.090 



2.1 X 10^ 

216 

0.41 

0.093 



2.3 X 10^ 

236 

0.28 

0.074 



2.6 X 10^ 

306 

0.29 

0.066 



6.3 X 10^ 

559 

0.13 

0.033 



9.7 X 10^ 

7.7 X 10® 

0.050 

8.7 X 10'® 



8.0 X 10"^ 

1.2 X 10® 

4.5 X 10-® 

1.3 X 10'® 



1.7 X 10^ 

925 

< 0.025 




+3.2 X 10^ 

2.3 X 10“* 

3.3 X 10^ 

1.6 X 10-+^ 

This work 


+2.0 X 10® 

3.4 X lO'* 

< 2.7 X 10^ 


This work 

141202A'' 

1.2 X 10® 

5.1 X 10“ 

< 6.1 X 10'® 


Pa^ani & Evans (20141 

WnOSA* 

2.9 X 10^ 

5.0 X 10® 

< 3.3 X 10'® 


^toling & Paee (20141 

141212A 

2.1 X 10“* 

2.9 X 10“ 

9.1 X 10'® 

3.5 X 10'® 

Evans et al. (2007a. 20091 

150101A 

1.3 X 10® 

2.3 X 10“ 

0.016 

3.7 X 10'® 

Evans et al. (2007a. 20091 


6.0 X 10“* 

1.8 X 10“ 

<0.010 



150101B 

+6.9 X 10® 

1.5 X 10“ 

0.011 

1.1 X 10'® 

Fong et al., in prep. 


+3.4 X 10® 

1.5 X 10“ 

1.9 X 10'® 

4.7 X 10^ 

Fong et al., in prep. 

150120A 

8.3 X 10® 

1.2 X 10“ 

<0.031 


Evans et al. (2007a. 20091 

150301A 

1.4 X 10^ 

4.0 X 10“ 

<0.051 


Evans et al. (2007a. 20091 


Note. — Upper limits correspond to 3cr. Unless otheiAvise stated, all data are taken with SwiftfXKI and X-ray flux densities ai‘e at 1 keV. 

° We employed a fiducial spectral index of = “1- 

^ We employed a fiducial count rate-to-unabsorbed flux conversion factor of 10“^^ erg cm“^ s“^ and spectral index [3x = —1- 
^ Reported flux is in the energy range of 0.7-10 keV. 

^ Late-time re-brightening in GRB 080503 light curve is observed in both optical and X-ray bands and is unlikely the afterglow tPerlev et al.l2009bl) . 
Chandra observation. 

XMM-Newton observation. 
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Table A2 

Short GRB Optical/near-IR afterglow catalog 


GRB 

St 

(hr) 

Telescope 

Instrument 

Lilter 

Exposure Time 
(S) 

(/iJy) 

(itJy) 

References 

050202 

12.6 

Mount John 

MOA 

R 

900 

< 42.89 


Castro-Tirado et al. 12005) 

050509B 

2.1 

WIYN 

OPTIC 

r 

600 

<0.75 


Bloom et a^2006) 

050709 

34.0 

Danish tel. 

DLOSC 

R 

7200 

2.34 

0.12 

Hjortli et 


59.1 

VLT 

PORS2 

V 

3600 

0.64 

0.07 

Covino et al. 120061 


59.3 

VLT 

PORS2 

R 

3000 

0.90 

0.05 



60.0 

Danish tel. 

DLOSC 

R 

10200 

1.17 

0.26 

Hiorth et al. (2005) 


104.7 

VLT 

LORSl 

V 

3600 

<0.36 


Covino et al.^20061 


134.4 

HST 

ACS 

L814W 

6360 

0.34 

0.006 

Lox et al. (20051 


235.2 

HST 

ACS 

L814W 

6360 

0.17 

0.008 



832.8 

HST 

ACS 

L814W 

6360 

<0.02 



050724A 

11.4 

VLT 

LORSl 

V 

480 

16.0 

0.45 

Malesani et al. (2007al 


11.6 

Magellan/Baade 

PANIC 

K 

1320 

45.6 

1.4 

Beraer et al. (20051 


11.8 

VLT 

LORSl 

R 

540 

17.8 

0.7 

Ma!esffli^^. (2nOTal 


11.8 

VLT 

LORSl 

I 

540 

18.37 

0.51 



12.0 

Swope 40-in 


I 

1800 

19.1 

0.2 

Berber et al. (20051 


14.2 

Swope 40-in 


I 

1800 

25.2 

0.9 



34.8 

VLT 

LORSl 

I 

540 

2.81 

0.33 

Malesani et al. (2007al 


34.9 

Magellan/Baade 

PANIC 

K 

1320 

<5.4 


Beraer et al. (20051 


35.0 

VLT 

LORSl 

R 

540 

3.35 

0.29 

Malesani et al. t2007al 


36.7 

Swope 40-in 


I 

2700 

< 9.1 


Beraer et al. (2(5(1^ 


83.1 

VLT 

LORSl 

I 

540 

0.33 

0.12 

Malesani et al. (2007al 

050813 

13.2 

CAHA 2.2-m 

CALOS 

I 

6000 

< 1.82 


Lerrero et al. (20071 


14.1 

CAHA 2.2-m 

CALOS 

R 

4140 

< 1.43 



050906 

21.4 

VLT 

PORS2 

R 

1800 

<0.15 


Levan et al^2008) 

051210 

19.4 

Magellan/Clay 

LDSS3 

r 

1200 

< 1.6 


Beraer et al. (20071 

051221A 

3.1 

Gemini-N 

GMOS 

r 


16.8 

1.3 

ioJeroere et^L^?no5) 


3.4 

Gemini-N 

GMOS 

r 


15.6 

1.2 



26.9 

Gemini-N 

GMOS 

r 


2.5 

0.19 



27.2 

Gemini-N 

GMOS 

i 


2.3 

0.54 



51.0 

Gemini-N 

GMOS 

i 


0.83 

0.32 



51.3 

Gemini-N 

GMOS 

z 


1.1 

0.47 



51.6 

Gemini-N 

GMOS 

r 


0.93 

0.11 



75.7 

Gemini-N 

GMOS 

r 


0.94 

0.28 



123.6 

Gemini-N 

GMOS 

r 


0.50 

0.11 



147.8 

Gemini-N 

GMOS 

r 


0.32 




060121 

0.37 

OSN 


/ 

120 

19.3 

4.4 

de Uaarte Postiao et al. 120061 


0.70 

OSN 


/ 

120 

10.0 

2.8 


2.0 

NOT 

ALFOSC 

R 


3.3 

0.70 

Levan et al. (2006) 


2.6 

NOT 

ALFOSC 

R 


6.1 

1.9 


2.7 

OSN 


I 

300 

10.8 

2.5 

de Uaarte Postiao et al. 120061 


2.8 

OSN 


R 

600 

3.7 

0.97 


3.9 

CAHA 2.2-m 


R 

1200 

1.1 

0.33 



4.4 

CAHA 2.2-m 


R 

1800 

4.6 

0.83 



5.4 

CAHA 2.2-m 


R 

1200 

1.2 

0.49 



5.6 

Bok 

90Prime 

R 


I.l 

0.22 

Levan et al. (2006) 


6.5 

Bok 

90Prime 

B 


<0.96 



7.3 

Bok 

90Prime 

R 


1.2 

0.18 



7.3 

CAHA 2.2-m 


R 

1800 

1.4 

0.37 



7.4 

WHT 


K 

750 

17.1 

1.4 

de Uaarte Postiao et al. 120061 


7.7 

OSN 


/ 

1500 

< 2.5 



11.3 

Bok 

90Prime 

R 


1.6 

0.41 

Levan et al. (2006) 


13.1 

WIYN 


R 


1.2 

0.24 


30.0 

OSN 


R 

10800 

< 0.8 


de Uaarte Postiao et al. 120061 


31.3 

WHT 


K 

1000 

6.3 

1.6 


33.5 

APO/ARC 

NIC-FPS 

K 

3600 

7.5 

0.65 

de Uaarte Postiao et al. t2006I 


33.6 

WIYN 


R 


0.41 

0.06 

Levan et al. (20061 


51.5 

CAHA 2.2-m 


R 

5400 

0.55 

0.14 

de U afflt^ostiTO e^T 120061 


127.8 

APO/ARC 

NIC-FPS 

K 

3600 

< 2.13 


060313 

0.045 

Swift 

UVOT 

V 

200 

55.0 

18.8 

Romina et al. 120061 


0.12 

Swift 

UVOT 

U 

50 

22.7 

13.0 


0.22 

Swift 

UVOT 

white 

50 

30.5 

11.6 



1.2 

Gemini-S 

GMOS 

r 

1800 

36.6 

6.6 

Berber et al. (20071 


1.3 

Swift 

UVOT 

B 

886 

18.7 

5.3 

Romina et al. 120061 


1.5 

Swift 

UVOT 

UVW2 

900 

6.4 

1.9 


1.7 

Swift 

UVOT 

V 

684 

33.7 

9.5 



2.9 

Swift 

UVOT 

UVM2 

900 

4.7 

2.3 



3.1 

Swift 

UVOT 

UVW\ 

900 

12.4 

2.4 



3.3 

Swift 

UVOT 

U 

708 

7.4 

1.8 



4.5 

Swift 

UVOT 

B 

886 

12.2 

5.8 



6.1 

Swift 

UVOT 

UVM2 

900 

3.4 

1.2 



6.3 

Swift 

UVOT 

UVWl 

900 

3.3 

1.5 



6.6 

Swift 

UVOT 

U 

716 

8.8 

2.1 
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Table A2 — Continued 


GRB 

St 

Telescope 

Instrument 

Filter 

Exposure Time 


<J 

References 


(hr) 




(S) 

(/iJy) 

(/rJy) 



9.6 

Swift 

UVOT 

UVWl 

900 

3.0 

1.4 



9.8 

Swift 

UVOT 

U 

658 

3.4 

1.9 



24.2 

Gemini-S 

GMOS 

r 

900 

3.7 

0.75 

Beraer et al. 120071 


48.5 

Gemini-S 

GMOS 

r 

1500 

1.3 

0.18 


240 

Gemini-S 

GMOS 

r 

1800 

<0.48 



060502B 

16.8 

Gemini-N 

GMOS 

R 

1500 

<0.62 


Price et al^2006I 

060801 

16.0 

Hale 

LEG 

r 

1500 

<0.83 


Beraer et al. t2007I 
D’Avanzo et al. (2UU9) 

061006 

14.9 

VLT 

FORSl 

I 

1800 

4.3 

0.20 


45.9 

VLT 

FORSl 

I 

1800 

2.5 

0.17 


60.0 

VLT 

FORSl 

I 

1260 

2.6 

0.22 



49.7 

VLT 

FORSl 

R 

645 

1.5 

0.17 


061201 

8.6 

VLT 

FORS2 

I 

3640 

2.9 

0.22 

Stratta et al. ('2007a^ 


9.0 

VLT 

FORS2 

R 

780 

1.9 

0.23 


33.1 

VLT 

FORS2 

I 

5200 

<0.91 




81.4 

VLT 

FORS2 

I 

2400 

<0.63 



061210 

2.1 

Gemini-N 

GMOS 

r 

5400 

< 1.44 


Beraer et al. 120071 

061217 

2.8 

Magellan/Clay 

LDSS3 

r 

300 

< 2.0 


070209 

0.29 

PROMPT 


V 

120 

< 81.66 


Johnson et al. f20071 

070406 

24.5 

NOT 

ALFOSC 

R 

1800 

< 1.26 


Malesani et al. (20(|)7h') 

070429B 

4.8 

Gemini-S 

GMOS 

R 

900 

<0.50 


Perley et a^20OTI 


27.1 

Blanco 

ISPI 

J 

1080 

<4.0 


Nvsewander et al. 12007') 

070707 

11.0 

VLT 

FORSl 

R 

1200 

2.1 

0.039 

Piranomonte et al. (2008) 


33.8 

VLT 

FORSl 

R 

1200 

1.0 

0.047 


37.4 

VLT 

FORSl 

R 

180 

0.82 

0.10 



59.3 

VLT 

FORSl 

R 

1200 

0.26 

0.020 



61.0 

VLT 

ISAAC 

J 

1800 

<0.60 




83.2 

VLT 

FORSl 

R 

3600 

0.078 

0.014 



108.0 

VLT 

FORSl 

R 

5100 

0.066 

0.015 


070714B 

0.21 

Liverpool 


r 

60 

37.2 

6.7 

Graham et al. (2009) 


0.27 

Liverpool 


r 

120 

43.1 

4.6 



0.30 

Liverpool 


i 

120 

50.7 

8.6 



0.40 

Liverpool 


r 

120 

37.9 

9.8 



0.44 

Liverpool 


r 

120 

38.3 

9.0 



23.6 

WHT 


R 

2400 

1.0 

0.22 



24.0 

TNG 

NICS 

J 


<9.1 


Covino et al. (2007) 

070724A 

2.3 

Gemini-N 

GMOS 

g 

360 

< 1.5 


Beraer et al. t2009I 


2.3 

Gemini-N 

GMOS 

i 

360 

1.1 

0.1 


2.8 

Gemini-N 

NIRI 

Ks 

900 

9.3 

1.5 



3.1 

Gemini-N 

NIRI 

J 

900 

3.4 

0.3 

Berjier et al. (299pL this work 


3.4 

Gemini-N 

NIRI 

H 

450 

7.8 

0.44 

Berger et al. (20091: this work 


3.7 

Gemini-N 

NIRI 

Ks 

900 

8.9 

1.5 

Beraer et al. 120091 

070729 

8.4 

ESO/MPG 

GROND 

r 


<0.48 


Nicuesa Guelbenzu et al. ('2012a) 


8.4 

ESO/MPG 

GROND 

J 


< 3.02 


070809 

11.2 

Keck I 

LRIS 

R 

640 

0.82 

0.17 

Perlev et al. I2007I 


35.1 

Keck I 

LRIS 

R 


0.36 

0.11 


11.21 

Keck I 

LRIS 

g 

880 

0.25 

0.051 



35.14 

Keck I 

LRIS 

g 


<0.19 



0708lOB 

23.4 

Keck I 

LRIS 

R 

630 

<0.21 


Kocevski et al. (20071 

071112B 

6.3 

Magellan/Clay 

LDSS3 

r 

1200 

< 2.0 


Rerper & Ghallis ('20()7) 


9.6 

ESO/MPG 

GROND 

J 


< 8.32 


Nicuesa Guelbenzu et al. (2012a) 

071227 

4.2 

ESO/MPG 

GROND 

J 


< 36.3 



7.0 

VLT 

FORS2 

R 

240 

1.6 

0.12 

D’Avanzo et al. (2009) 


80.9 

VLT 

FORS2 

R 

540 

0.86 

0.15 



080121 

55.2 

Swift 

UVOT 

white 

2015 

<4.57 


Troia et al. (20081 

080123 

0.04 

Swift 

UVOT 

white 

100 

<22.91 


Ukwatta et al. (20081 

080426" 

7.5 

CAHA 2.2-m 


I 

3300 

<7.06 


de Ilaarte Posfiao et al. 120081 

080503 

0.037 

Swift 

UVOT 

white 

98 

< 14.2 


Perlev et al. (200961 


1.0 

Gemini-N 

GMOS 

r 

180 

<0.20 




1.2 

Gemini-N 

GMOS 

g 

900 

0.089 

0.018 



1.2 

Keck I 

LRIS 

B 

300 

<0.21 




1.3 

Keck I 

LRIS 

R 

630 

<0.21 




1.5 

Gemini-N 

GMOS 

r 

800 

<0.081 




1.8 

Gemini-N 

GMOS 

i 

800 

< 0.078 




2.2 

Gemini-N 

GMOS 

z 

800 

<0.16 




2.4 

Gemini-N 

GMOS 

g 

360 

<0.65 




*26.0 

Gemini-N 

GMOS 

r 

1800 

0.27 

0.037 



*47.4 

Gemini-N 

GMOS 

r 

1620 

0.23 

0.038 



*50.2 

Gemini-N 

GMOS 

g 

720 

0.12 

0.024 



*74.0 

Gemini-N 

GMOS 

r 

2700 

0.19 

0.046 



*97.1 

Gemini-N 

GMOS 

r 

2880 

0.13 

0.025 



125 

Gemini-N 

NIRI 

Ks 

2760 

<0.70 
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GRB 

St 

(hr) 

Telescope 

Instrument 

Filter 

Exposure Time 
(S) 

(/iJy) 

(/rJy) 

References 


*128.6 

HST 

WFPC2 

F606W 

0.067 

0.011 



080702A 

12.06 

Loiano 


R 

1800 

< 44.92 


Greco et al. (20081 

080905A 

8.5 

NOT 

ALFOSC 

R 

1800 

0.72 

0.39 

Rowlinson et al. f2010bl 


14.3 

VLT 

FORS2 

R 

2400 

0.59 

0.20 



17.5 

ESO/MPG 

GROND 

r 

660 

< 3.7 


Nicuesa Guelbenzu et al. t2012a) 


17.5 

ESO/MPG 

GROND 

j 

660 

<27.5 




36.0 

VLT 

FORS2 

R 

2400 

<0.30 


Rowlinson et al. 1201061 

080919“ 

0.19 

ESO/MPG 

GROND 

J 


<133 


Nicuesa Guelbenzu et al. (2012al 

081024A 

1.9 

Faulkes North 


R 

600 

< 93.93 


Melandri et al. (20081 

081024B 

30 

P200 

LFC 


R 

< 1.36 


Genko & Kasliwal ('20081 

081226A 

0.37 

ESO/MPG 

GROND 

r 


1.9 

0.42 

Nicuesa Guelbenzu et al. f2012al 


1.1 

ESO/MPG 

GROND 

g 


0.38 

0.12 



1.1 

ESO/MPG 

GROND 

r 


0.63 

0.16 



1.1 

ESO/MPG 

GROND 

i 


0.81 

0.31 



1.1 

ESO/MPG 

GROND 

z 


1.4 

0.35 



1.1 

ESO/MPG 

GROND 

J 


<7.7 



081226B 

2.7 

Swift 

UVOT 

B 

1414 

< 14.45 


Evans & Floversten 120081 

090305 

0.47 

ESO/MPG 

GROND 

r 


3.0 

0.084 

Nicuesa Guelbenzu et al. t2012al 


0.56 

ESO/MPG 

GROND 

r 


3.1 

0.27 



0.56 

ESO/MPG 

GROND 

i 


3.3 

0.59 



0.60 

ESO/MPG 

GROND 

r 


2.9 

0.082 



0.61 

Gemini-S 

GMOS 

r 

900 

1.8 

0.034 

Tunnicliffe et al. (20131 


0.71 

ESO/MPG 

GROND 

g 


2.1 

0.40 

Nicuesa (jueltenzu et al. ('5oi2a) 


0.71 

ESO/MPG 

GROND 

r 


2.8 

0.35 



0.79 

ESO/MPG 

GROND 

g 


1.9 

0.09 



0.92 

ESO/MPG 

GROND 

g 


2.1 

0.31 



0.92 

ESO/MPG 

GROND 

r 


2.1 

0.32 



0.92 

ESO/MPG 

GROND 

i 


2.5 

0.49 



0.93 

ESO/MPG 

GROND 

g 


1.8 

0.085 



1.1 

ESO/MPG 

GROND 

g 


1.7 

0.11 



1.1 

ESO/MPG 

GROND 

r 


1.9 

0.24 



1.2 

ESO/MPG 

GROND 

r 


2.0 

0.13 



1.3 

Gemini-S 

GMOS 

i 

900 

1.5 

0.041 

Tunnicliffe et al. (20131 


1.3 

ESO/MPG 

GROND 

g 


1.6 

0.14 

Nicuesa Guelnenzu et al. (z012al 


1.3 

ESO/MPG 

GROND 

i 


1.9 

0.24 



1.3 

ESO/MPG 

GROND 

r 


1.5 

0.16 



1.4 

ESO/MPG 

GROND 

r 


1.7 

0.082 



1.5 

ESO/MPG 

GROND 

r 


1.7 

0.34 



1.5 

ESO/MPG 

GROND 

i 


1.9 

0.22 



1.6 

ESO/MPG 

GROND 

r 


1.7 

0.078 



1.6 

Gemini-S 

GMOS 

r 

900 

1.4 

0.026 

Tunnicliffe. et al. (2(11^1 


1.6 

ESO/MPG 

GROND 

r 


1.8 

0.19 

Nicuesa Guelbenzu et al. ('2012al 


1.7 

ESO/MPG 

GROND 

r 


1.6 

0.073 



1.7 

ESO/MPG 

GROND 

r 


1.5 

0.12 



1.8 

ESO/MPG 

GROND 

r 


1.4 

0.051 



1.8 

ESO/MPG 

GROND 

g 


1.0 

0.13 



1.8 

ESO/MPG 

GROND 

r 


1.5 

0.12 



2.0 

ESO/MPG 

GROND 

r 


1.6 

0.12 



2.1 

ESO/MPG 

GROND 

r 


1.2 

0.35 



2.1 

Gemini-S 

GMOS 

r 

900 

0.95 

0.027 

Tunnicliffe et al. (20131 


2.1 

ESO/MPG 

GROND 

r 


1.1 

0.040 

Nicuesa Guelbenzu et al. ('2012al 


21.7 

Gemini-S 

GMOS 

r 

1500 

<0.19 


Tunnicliffe et al. (20131 

090426 

0.29 

TNT 


R 

60 

101 

7.8 

Xinet al. (20111 


0.31 

TNT 


R 

60 

95.1 

lO.l 



0.33 

TNT 


R 

60 

89.2 

9.5 



0.36 

TNT 


R 

60 

77.0 

8.2 



0.38 

TNT 


R 

60 

74.9 

8.8 



0.40 

TNT 


R 

60 

69.6 

8.1 



0.42 

TNT 


R 

60 

73.5 

8.6 



0.53 

TNT 


V 

600 

49.0 

5.7 



0.71 

TNT 


V 

600 

36.5 

4.3 



0.88 

TNT 


V 

600 

34.5 

4.0 



l.I 

TNT 


V 

600 

36.2 

4.2 



1.2 

TNT 


V 

600 

32.4 

3.8 



1.4 

TNT 


V 

600 

31.8 

4.4 



1.6 

TNT 


V 

600 

28.7 

4.0 



2.0 

TNT 


V 

600 

24.6 

3.9 



2.3 

TNT 


V 

600 

20.1 

3.4 



2.6 

TNT 


V 

600 

20.8 

3.5 



3.0 

TNT 


V 

600 

19.7 

3.6 



7.3 

TLS 


I 


15.8 

2.8 

Nicuesa Guelbenzu et al. (20111 


7.5 

TLS 


R 


13.5 

1.7 



7.5 

TLS 


R 


13.4 

1.8 
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Table A2 — Continued 


GRB 

St 

(hr) 

Telescope 

Instrument 

Filter 

Exposure Time 
(S) 

(/iJy) 

(/rJy) 

References 



7.7 

TLS 


R 


11.8 

2.3 




7.8 

TLS 


R 


12.7 

1.7 




12.9 

ESO/MPG 

GROND 

g 

3.1 xlO^ 

5.8 

0.16 




12.9 

ESO/MPG 

GROND 

r 

3.1 xlO^ 

7.6 

0.21 




12.9 

ESO/MPG 

GROND 

i 

3.1 xlO^ 

8.2 

0.38 




12.9 

ESO/MPG 

GROND 

z 

3.1 xlO^* 

9.8 

0.65 




12.9 

ESO/MPG 

GROND 

j 

3.1 xlO^ 

15.2 

0.28 




12.9 

ESO/MPG 

GROND 

H 

3.1 xlO^* 

15.7 

3.2 




12.9 

ESO/MPG 

GROND 

K 

3.1 xlO^ 

< 33.3 





14.6 

ESO/MPG 

GROND 

8 

7.3 X 10^ 

4.7 

0.087 




14.6 

ESO/MPG 

GROND 

r 

7.3 X 10^ 

6.0 

0.11 




14.6 

ESO/MPG 

GROND 

i 

7.3 xl03 

7.3 

0.20 




14.6 

ESO/MPG 

GROND 

z 

7.3 X 10^ 

8.3 

0.31 




14.6 

ESO/MPG 

GROND 

j 

7.3 xl03 

10.2 

0.89 




14.6 

ESO/MPG 

GROND 

H 

7.3 X 10^ 

11.8 

1.4 




14.6 

ESO/MPG 

GROND 

K 

7.3 X 10^ 

<21.0 





16.3 

LOAO 


R 

1.8 xlO^ 

< 6.6 


Xinetal. (20111 

090510'' 

6.2 

ESO/MPG 

GROND 

r 


5.9 

2.5 

Nicuesa Guelbenzu et al. 12012b) 


6.3 



r 


5.5 

2.3 




6.3 



r 


4.6 

1.6 




6.4 



r 


2.6 

1.8 




6.4 



r 


3.0 

1.8 




6.5 



r 


2.6 

1.6 




6.7 



r 


2.3 

0.71 




6.8 



r 


3.8 

0.77 




6.9 



r 


2.7 

0.89 




7.1 



r 


3.4 

0.79 




7.2 



r 


2.9 

0.72 




7.3 



r 


2.2 

0.75 




7.4 



r 


3.1 

0.77 




7.6 



r 


2.3 

0.59 




7.7 



r 


1.9 

0.40 




7.8 



r 


1.9 

0.50 




7.9 



r 


2.2 

0.55 




8.1 



r 


1.8 

0.59 




8.2 



r 


2.0 

0.52 




8.3 



r 


2.0 

0.47 




8.4 



r 


2.1 

0.46 




8.6 



r 


2.3 

0.48 




8.6 



J 


<4.8 





8.7 



r 


2.0 

0.47 




8.8 



r 


1.9 

0.44 




8.9 



r 


1.5 

0.47 




9.1 



r 


1.7 

0.43 




9.2 



r 


1.3 

0.41 




9.3 



r 


1.5 

0.41 




9.5 



r 


1.5 

0.40 




9.7 



r 


0.88 

0.47 




9.8 



r 


1.2 

0.41 



090515 

1.4 

WlYN 

WHIRC 

J 

2.4 X 103 

<39.4 


Updike et al. (2009) 


1.7 

Gemini-N 

GMOS 

r 

1.8 xlO^ 

0.11 

0.013 

Rowlinson et al. ('2010a) 


2.0 

MMT 

MMIRS 

K 

810 

< 12.0 


McLeod & Williams (2009) 


25.0 

Gemini-N 

GMOS 

r 

1.8 xlO^ 

0.092 

0.033 

Rowlinson et al. (2010al 


4.4 X 10^ 

Gemini-N 

GMOS 

r 

2.8 xlO^ 

< 0.042 




090607“ 

0.52 

Faulkes North 


R 

720 

< 5.70 


(jriiidor/i et al. (2()()9) 

09062 IB 

0.75 

RTT150 

TFOSC 

Rc 

900 

< 1.91 


Galeev et al. (2009) 

090916 

2.7 

PROMPT 


R 

480 

<71.4 


Haislin et al. (20091 

091109B 

6.0 

VET 

FORS2 

R 

2400 

0.67 

0.092 

Tunnicliffe et al. (2013) 


7.2 

VET 

HAWK-1 

K 

1320 

<0.85 





8.3 

VET 

HAWK-1 

J 

1320 

< 2.5 





10.4 

VET 

FORS2 

R 

1200 

0.51 

0.11 




31.7 

VET 

FORS2 

R 

2400 

<0.17 




091117 

31.5 

ESO/MPG 

GROND 

J 


<7.59 


Nicuesa Guelbenzu et al. (2012a) 


31.7 

Magellan 

IMACS 

r 


<0.98 


Beraer & Mulchaev 120091 

100117A 

3.9 

Magellan 

IMACS 

R 

1200 

<0.97 


Fona et al. (20111 


4.3 

ESO/MPG 

GROND 

r 


<0.69 


Nicuesa Guelbenzu et a 

(N 

O 

d 

1 -1 


8.3 

Gemini-N 

GMOS 

r 

2700 

0.24 

0.05 

Fona et al. (201 

D 

100206A 

11.7 

ESO/MPG 

GROND 

i 


< 1.8 


Nicuesa Guelbenzu et al. (2012a) 


11.7 

ESO/MPG 

GROND 

J 


< 9.7 





15.7 

Gemini-N 

GMOS 

i 

1200 

<0.18 


Perley et al. (20121 

100625A 

12.2 

ESO/MPG 

GROND 

g 

fa 3600 

< 1.32 


Nicuesa Guelbenzu et al. (2012a) 
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GRB 

St 

(hr) 

Telescope 

Instrument 

Eilter 

Exposure Time 
(S) 

(/iJy) 

(itJy) 

References 


33.9 

Magellan 

PANIC 

J 

2100 

< 1.0 


Fonj; et al^^ii^' 

100628A 

1.1 

Gemini-N 

GMOS 

i 

1200 

<0.85 


This work 


17.8 

Magellan 

PANIC 

j 

1620 

<5.43 



100702A" 

0.16 

ESO/MPG 

GROND 

j 


< 27.29 


Nicuesa Guelbenzu et al. t2012a) 


1.7 

ESO/MPG 

GROND 

J 


r 

< 1.96 


101219A 

0.96 

Gemini-S 

GMOS 

i 

1620 

<0.40 


Fong et al. Q0131 


1.20 

Magellan 

FourStar 

J 

1500 

< 1.36 



110112A 

15.4 

WHT 

ACAM 

i 

600 

2.8 

0.75 

Fong et al. (201.31 


4.0 X10^ 

Magellan 

LDSS3 

i 

1200 

<0.49 



110112B 

17.0 

Swift 

UVOT 

white 

3600 

< 5.01 


Littleiohns et al. (201 1) 

110420B 

10.5 

Magellan 

IMACS 

r 

1530 

< 1.4 


This work 

111020A“ 

17.8 

Magellan/Clay 

LDSS3 

r 

1080 

<0.83 


Fong et al. (2012bl 


17.8 

Gemini-S 

GMOS 

i 

1620 

<0.63 




18.5 

VET 

HAWK-I 

J 

2640 

<0.57 


Tunnicliffe et al. (2013) 

111117A 

8.2 

GTC 

OSIRIS 

r 

1200 

<0.40 


Sakamoto et al. (2013) 


13.7 

Gemini-S 

GMOS 

r 

1200 

<0.23 


Marautti et ^^(21?^ 


120229A 

9.6 

Magellan/Clay 

LDSS3 

r 

540 

<4.1 

This work 

120305A 

0.23 

Liverpool 


r 


< 25.85 

Vireili et al. 120121 

120521A“ 

18.6 

ESO/MPG 

GROND 

r 

3000 

<0.88 

Rossi et al. 120121 


18.6 

ESO/MPG 

GROND 

J 

2400 

< 20.9 


120804A 

5.5 

Gemini-N 

GMOS 

i 

1980 

0.17 0.04 

Beraer et al. 120131 

120817B 

24.0 

LCO/duPont 

WFCCD 

R 

900 

< 1.5 

Fona et al. 12012al 

121226A 

11.1 

NOT 

ALFOSC 

R 

1800 

< 2.1 

Xuetal. 120121 


130313A“ 

13.4 

TNG 


r 

1800 

<0.44 


D’Avanzn et al. 120131 

130515A 

0.63 

Gemini-S 

GMOS 

r 

240 

< 1.4 


Cenko & Cucchiara (2013) 

130603B 

0.19 

Swift 

UVOT 

V 

180 

< 199.5 


de Uaarte Postiao et al. (20131 


2.1 

Swift 

UVOT 

V 

5110 

<54.0 




2.4 

Swift 

UVOT 

B 

7700 

17.2 

4.5 



5.9 

NOT 

MOSCA 

r 

1800 

12.6 

0.23 



6.1 

WHT 

ACAM 

z 

900 

25.4 

1.4 



6.6 

WHT 

ACAM 

i 

900 

16.4 

0.88 



6.7 

CAHA 

DLR-MKIII 

V 

1800 

8.3 

0.73 



7.0 

GTC 

OSIRIS 

r 

30 

11.0 

0.20 



7.1 

WHT 

ACAM 

g 

900 

6.3 

0.34 



7.9 

Gemini-S 

GMOS 

g 

1440 

5.3 

0.19 

Cucchiara et al. (2013) 


8.2 

Magellan/Baade 

IMACS 

r 

600 

8.6 

0.14 

Beraer et al. (20131 


9.0 

Gemini-S 

GMOS 

i 

1440 

12.3 

1.2 

Cucchiara et al. (2013) 


14.4 

Gemini-N 

GMOS 

z 

500 

6.5 

0.18 

de Uaarte Postiao et al. (201.31 


14.5 

UKIRT 

WFCAM 

K 

684 

13.7 

1.3 



14.6 

Gemini-N 

GMOS 

i 

500 

4.5 

0.12 



14.7 

UKIRT 

WFCAM 

J 

500 

9.3 

1.3 



14.8 

Gemini-N 

GMOS 

r 

500 

2.9 

0.081 



15.0 

Gemini-N 

GMOS 

g 

500 

1.6 

0.06 



31.2 

Gemini-S 

GMOS 

r 

540 

<0.30 


Cucchiara et al. (2013) 


31.2 

Gemini-S 

GMOS 

i 

540 

<0.58 




32.2 

Magellan/Baade 

IMACS 

r 

1200 

<0.46 


Berger et al. (20131 


38.2 

Gemini-N 

GMOS 

g 

600 

<0.19 


de Uaarte Postiao e^n'20131 


38.4 

Gemini-N 

GMOS 

r 

600 

0.21 

0.05 



38.6 

Gemini-N 

GMOS 

i 

600 

<0.48 




38.6 

UKIRT 

WFCAM 

J 

1400 

<3.6 



130716A 

19.5 

Gemini-N 

GMOS 

r 

900 

<0.35 


This work 

130822A 

21.1 

Gemini-N 

GMOS 

i 

600 

<0.58 


Cenko et al. (2013a) 

130912A 

0.27 

ESO/MPG 

GROND 

r 


10.5 

2.1 

Tan^a et al. (2013) 


0.38 

P60 


r 


13.0 

2.6 

Cenko et al. (2013b) 


0.88 

P60 


r 


9.7 

2.5 



20.0 

WHT 

ACAM 

g 


<0.87 


Tanvir et al. (20131 


24.8 

Johnson 

RATIR 

r 

19200 

< 1.8 


Butler et al. (20131 


24.8 

Johnson 

RATIR 

i 

19200 

< 1.7 




24.8 

Johnson 

RATIR 

z 

8060 

<3.7 




24.8 

Johnson 

RATIR 

Y 

8060 

< 5.4 




24.8 

Johnson 

RATIR 

J 

8060 

< 5.3 




24.8 

Johnson 

RATIR 

H 

8060 

< 7.9 



131004A 

0.43 

NOT 

ALFOSC 

R 


44.3 

4.3 

Xu et al. (20131 


1.1 

TNG 

DOLORES 

R 

60 

25.5 

2.5 

MaJesmi^ a^55l.31 


1.4 

Liverpool 


i 

300 

47.4 

4.6 

Kopac et al. (2013) 


1.9 

Liverpool 


i 

300 

32.8 

3.2 



2.8 

ESO/MPG 

GROND 

g 

460 

6.0 

0.58 



2.8 

ESO/MPG 

GROND 

r 

460 

9.6 

0.93 



2.8 

ESO/MPG 

GROND 

i 

460 

13.1 

1.3 
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Table A2 — Continued 


GRB 

St 

(hr) 

Telescope 

Instrument 

Filter 

Exposure Time 
(S) 

(/iJy) 

(/tJy) 

References 


2.8 

ESO/MPG 

GROND 

z 

460 

14.0 

1.4 



2.8 

ESO/MPG 

GROND 

J 

480 

36.5 

7.4 



2.8 

ESO/MPG 

GROND 

H 

480 

34.0 

6.9 



2.8 

ESO/MPG 

GROND 

K 

480 

<47.3 




13 

Johnson 

RATIR 

r 

7812 

< 1.6 


Littleiohns et al. (2013) 


13 

Johnson 

RATIR 

i 

7812 

< 1.9 




13 

Johnson 

RATIR 

z 

3384 

< 2.9 




13 

Johnson 

RATIR 

Y 

3384 

<4.8 




13 

Johnson 

RATIR 

J 

3384 

<3.6 




13 

Johnson 

RATIR 

H 

3384 

<3.6 



131125A“ 

11.8 

iPTF 


R 


<24.9 


Sinper et al. (2()1 Jla) 

131126A'' 

19.5 

iPTF 


R 


< 107.8 


Sinser et al. 12013b) 

140129B 

0.032 

MASTER II 



20 

7590 


Poleshchuk et al. (2014) 


0.053 

MASTER II 


- 

30 

3020 




0.079 

MASTER II 


. 

50 

1910 




0.11 

MASTER II 


- 

70 

1450 




0.15 

MASTER II 


- 

100 

692 




0.20 

MASTER II 


- 

130 

525 




0.26 

MASTER II 


- 

170 

363 




0.33 

MASTER II 


- 

180 

275 




0.43 

MASTER II 


- 

360 

191 



140320A 

32.9 

Johnson 

RATIR 

r 

5040 

< 10.7 


Littleiohns et al. (2014) 


32.9 

Johnson 

RATIR 

i 

5040 

< 10.4 




32.9 

Johnson 

RATIR 

z 

2016 

<24.3 




32.9 

Johnson 

RATIR 

Y 

2016 

<20.2 




32.9 

Johnson 

RATIR 

J 

2016 

< 25.3 




32.9 

Johnson 

RATIR 

H 

2016 

< 37.3 



140402A" 

29.0 

Magellan/Baade 

IMACS 

r 

900 

<0.36 


This work 

140414A 

15.2 

TNG 

DOLORES 

r 

1800 

< 5.9 


D’Avanzo et al. (2014a) 

140516A 

1.6 

TNG 

DOLORES 

r 

720 

<2.38 


D’Avanzo et al. (2014b) 


1.8 

NOT 

ALFOSC 

R 

3900 

<0.50 


Gorosabel et al. (2014b) 


11.6 

Subaru 

IRCS+A0188 

K' 

5880 

<0.84 


Minowa et al. (2014) 


12.4 

Gemini-N 

GMOS 

i 

1800 

<0.14 


This wo^^^" 


13.0 

Johnson 

RATIR 

J 

3200 

< 5.8 


Butler et al. (20141 

140606A 

7.3 

AAO 


R 

120 

<6.46 


Voinova et al^20141 

140619B” 

35.5 

Magellan/Baade 

Fourstar 

J 

770 

<2.54 


This work 

140622A 

0.07 

ESO/MPG 

GROND 

r 


<0.80 


Tan^a et al. (20141 


0.07 

ESO/MPG 

GROND 

J 


< 11.4 



140903A” 

12.2 

DCT 

LMI 

r 


7.6 

3.1 

Cenko et a^20H) 

140930B 

3.0 

WHT 

ACAM 

r 


7.4 


Tanvir et al. (20141 


7.3 

ESO/MPG 

GROND 

r 


2.0 

0.41 

Graham et al. (2014) 


7.8 

Magellan/Baade 

IMACS 

i 

900 

2.2 

0.30 

This work 


7.8 

MMT 

MMTCam 

r 

1200 

0.97 

0.20 

This work 


13.1 

Keck 

MOSFIRE 

Ks 

315 

3.1 

0.97 

Perlev & Jencson (2014) 


13.3 

Keck 

MOSFIRE 

J 

396 

2.1 

0.43 


141205A 

11.4 

Nanshan 


R 


< 33.2 


Xu et al. (20141 

141212A 

5.4 

CAHA 1.23-m 


I 


< 11.1 


Gorosabel et al. (2014a) 

150101B 

39.8 

Magellan/Baade 

IMACS 

r 

1200 

2.26 

0.34 

Fong et al., in prep. 


63.6 

Magellan/Baade 

IMACS 

r 

1200 

1.4 

0.38 



231 

TNG 

NICS 

J 


<7.59 


D’Avanzo et al. (20151 


257 

Gemini-South 

GMOS 

r 

1710 

<0.76 


Fong et al., in prep. 


349 

TNG 

NICS 

J 


< 3.98 


D’Avanzo et al. (20151 


360 

VET 

HAWK-I 

H 


< 1.45 


van der Horst et al. (i6l51 


694 

Magellan/Baade 

IMACS 

r 

1200 

<0.53 


Fong et al., in prep. 

150120A 

2.0 

Gemini-North 

GMOS 

r 

1800 

<0.36 


This work 


Note. — Upper limits correspon d to 3cr. All optical and near-IR fluxes are corrected for Galactic extinction in the directions of each bursts tSchle 2 el et alJ 
ll998l : ISchlaflv & Finkbein^l201 ID . Instruments, exposure times, and flux density uncertainties are provided whenever available. 

° Optical observing constraint, due to delayed precise localization, sightline with high Galactic extinction, crowded field, or proximity to a bright star which 
contaminates the GRB position. 

^ Optical re-brightening observed in GRB 080503 is unli kely the afterglow jPerley et al.|20 09bi) . 

^ Full simultaneous giz-band light curves are available in INicuesa Guelbenzu et ain^O^H) . 
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Table A3 

Short GRB Radio afterglow catalog 


GRB 

St 

Facility Mean Frequency 

Fi, a 

References 


(dy) 

(GHz) 

(fiJy) (fiJy) 



050202 

0.52 

VLA 

4.86 

< 120 


Frail & Soderbers (2005) 

050509B 

0.49 

WSRT 

4.86 

< 66 


van der Horst et al. (2005) 

050709 

1.6 

VLA 

8.46 

< 115 


Fox et al. (20051 


2.5 

VLA 

8.46 

< 114 




4.5 

VLA 

8.46 

<74 




7.5 

VLA 

8.46 

<40 



050724A 

0.57 

VLA 

8.46 

173 

30 

Berser et al. (2005) 


1.68 

VLA 

8.46 

465 

29 



9.10 

VLA 

8.46 

<259 


Panaitescu 120061 

050813 

1.64 

VLA 

8.46 

<55 


Cameron & Frail 12005al 

050906 

3.91 

VLA 

8.46 

<92 


Cameron & Frail (2005bl 

050925 

0.40 

WSRT 

4.9 

<72 


van der Horst (20051 

051105A 

0.54 

VLA 

8.5 

<51 


Frail & Cameron (20051 

051221A 

0.91 

VLA 

8.46 

155 

30 

Soderbers et al. (20061 


1.94 

VLA 

8.46 

<72 




060313 

2.12 

VLA 

8.46 

< 110 

Soderbers & Frail 120061 

060801 

0.49 

VLA 

8.46 

< 105 

Soderbers et al. 120061 


1.2 

WSRT 

4.9 

<72 

van der Horst 120061 


5.2 

WSRT 

4.9 

<81 

061210 

1.9 

VLA 

8.46 

< 102 

Chandra & Frail (20061 


070429B 

0.59 

VLA 

8.46 

<300 

Chandra & Frail (2007bl 

070714B 

15.6 

VLA 

8.46 

< 135 

Chandra & Frail (2007dl 

070724A 

1.06 

VLA 

8.46 

<255 

Chandra & Frail (2007al 

070729 

9.49 

VLA 

8.46 

<255 

Chandra & Frail (2007el 

070923 

4.98 

VLA 

8.46 

< 135 

Chandra & Frail (2007cl 

071112B 

1.69 

ATCA 

8.7 

< 141 

Wierinsa et al. (20071 


080503 

3.05 

VLA 

8.46 

<54 

Frail & Chandra (2008) 

080702A 

0.74 

VLA 

8.46 

< 156 

Chandra & Frail (2008b) 

081024A 

1.03 

VLA 

8.46 

< 156 

Chandra & Frail (2008a) 

081024B 

4.1 

VLA 

8.46 

< 114 

Chandra et al. (2008) 

081226A 

58.9 

ATCA 

4.9 

<540 

Moin et al. 12009al 

081226B 

58.4 

ATCA 

4.9 

<588 

Moin et al. 12009b) 

090417A 

0.37 

VLA 

8.46 

<72 

Chandra & Frail 120091 

090510 

1.98 

VLA 

8.46 

<84 

Frail & Chandra 120091 

090515 

0.87 

VLA 

8.46 

<60 

Berser & Fons 12009al 

09062IB 

0.61 

VLA 

8.46 

<54 

Berser & Fons 12009cl 

090715A 

1.25 

VLA 

8.46 

<64 

Berser & Fons (2009d) 

091117 

2.33 

VLA 

8.46 

< 120 

Berser & Fons (2009bl 

100625A 

0.83 

VLA 

4.90 

< 192 

This work 

100628A 

0.77 

VLA 

5.8 

<291 

This work 

101224A 

3.72 

VLA 

5.8 

<56 

This work 

110112A 

1.87 

VLA 

5.0 

<75 

This work 

110112B 

0.85 

VLA 

5.8 

<51 

This work 


1.80 

VLA 

5.8 

<66 



2.92 

VLA 

5.8 

<48 



12.10 

VLA 

5.8 

<36 


110420B 

0.65 

VLA 

4.90 

<77 

This work 


7.65 

VLA 

5.8 

< 116 


111020A 

0.67 

VLA 

5.8 

<39 

Fons et al. 12012bl 

111117A 

0.49 

VLA 

5.8 

< 18 

Marsutti et al. 120121 

111121A 

0.84 

VLA 

5.8 

< 100 

This work 

120229A 

0.38 

VLA 

5.8 

<54 

This work 


3.32 

VLA 

5.8 

< 66 


120305A 

0.16 

VLA 

5.8 

< 18 

This work 


2.25 

VLA 

5.8 

< 18 


120521A 

1.08 

ATCA 

34 

<95 

Hancock et al. (2012al 

120804A 

0.9 

VLA 

5.8 

<20 

Berser et al. 12013) 


2.2 

ATCA 

34 

<200 

Hancock et al. (2012b) 


4.2 

ATCA 

34 

< 120 


121226A 

1.75 

VLA 

5.8 

<30 

This work 


4.76 

VLA 

5.8 

<45 


130313A 

0.77 

VLA 

5.8 

<60 

This work 

130603B 

0.37 

VLA 

4.9 

125 14.4 

Fons et al. (2014) 
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Table A3 — Continued 


GRB 

St 

(dy) 

Facility 

Mean Frequency 
(GHz) 

F, 

(^Jy) 

a 

(^Jy) 

References 


0.37 

VLA 

6.7 

119 

9.1 



1.43 

VLA 

4.9 

<57 




1.43 

VLA 

6.7 

64.9 

15.2 



1.44 

VLA 

21.8 

<50 




4.32 

VLA 

4.9 

<51 




4.32 

VLA 

6.7 

<26 




84.31 

VLA 

4.9 

<69 




84.31 

VLA 

6.7 

<34 



130716A 

2.51 

VLA 

5.8 

<33 


This work 

130822A 

15.9 

VLA 

5.8 

<30 


This work 

130912A 

0.91 

VLA 

5.8 

<29 


This work 


3.10 

VLA 

5.8 

<42 



131004A 

0.14 

CARMA 

93 

< 300 


This work 

131224A 

0.20 

VLA 

6.0 

<33 


This work 


2.27 

VLA 

6.0 

<33 



140516A 

0.62 

VLA 

6.0 

<24 


This work 


4.67 

VLA 

6.0 

<21 



140619B 

1.51 

VLA 

6.0 

<36 


This work 

140622A 

3.09 

VLA 

6.0 

<36 


This work 

140903A 

0.40 

VLA 

6.0 

110 

9.5 

This work 


2.45 

VLA 

6.0 

187 

8.7 



2.45 

VLA 

9.8 

151 

8.9 



4.70 

VLA 

6.0 

128 

15.1 



7.42 

GMRT 

1.39 

102 

33 

Navana & Chandra (2014) 


9.24 

VLA 

6.0 

81.9 

14.7 

This work 


9.24 

VLA 

9.8 

<69 




18.24 

VLA 

6.0 

< 120 



140930B 

0.29 

VLA 

9.8 

<25.5 


This work 


3.27 

VLA 

9.8 

< 24.4 



141212A 

0.45 

VLA 

6.0 

<25.2 


This work 


3.76 

VLA 

6.0 

27.0 

8.1 



7.72 

VLA 

6.0 

21.9 

8.3 



7.72 

VLA 

9.8 

21.3 

6.4 



76.49 

VLA 

6.0 

< 19.5 



150101A 

1.43 

VLA 

9.8 

<48.6 


This work 

150101B 

5.73 

VLA 

9.8 

< 375 


Fong et ah, in prep. 

150120A 

0.90 

VLA 

9.8 

<29.6 


This work 


Note, 


For bursts with multiple contemporaneous upper limits, we only display the most constraining limits. Upper limits are 3cr. 









